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I. INFRODUCTION 


A. CFAR BACKGROUND 

The primary task of any radar system is to detect all objects in some area of obser- 
vation and to estimate their position. One current trend in carrying out this task is au- 
tomatic detection and tracking systems. This trend is driven by the demands of users 
who place a premium on data acquisition systems which provide rapid, highly accurate 
information. The basic precept of an automatic detection radar is the elimination of the 
human operator who was solely responsible for target detection. Early radar systems 
were designed to route all incoming information directly to the users video display. 
Clutter, noise and target amplitude variations were all displayed simultaneously. Target 
detection was therefore delegated to a highly trained operator who usually distinguished 
targets from interfering background noise, clutter and possibly electronic jamming. The 
requirement to replace the human operator involved factors such as operator fatigue, 
Saturation and reaction time. Studies have shown that an operator can simultaneously 
track only a few targets accurately for any extended period of time. As the weakest link 
in the target detection and tracking problem, the human operator has been replaced by 
advanced digital signal processing technologies that satisfy the intense requirements of 
todays radar users. 

Automatic target detection and tracking would be a simple task if the echoing object 
was always located in the clear. In this case, the signal received could simply be com- 
pared with some fixed threshold and targets declared whenever the signal exceeds this 
threshold. In actual application however, the target generally appears before a back- 
ground filled with complicated clutter tvpes and various sources of noise interference 
including jamming energies. It is clear that an automatic radar must have some means 
of false target rejection that was previously delegated to the skilled operator. The signal 
processing system that completes this task is known as the Constant False Alarm Rate 
(CFAR) processor. Simply put, the CFAR processor is an algorithm used by automatic 
detection radars to control the rate at which target detections are falsely declared. The 
ideal CFAR detector would be one which maximizes the probability of detection of a 
target when it does appear and minimizes the probability of false alarm caused by noise 
and clutter when no target is present. Unfortunately both problems cannot be optimized 


simultaneously. The best the CFAR algorithm can do is keep the false alarm rate at 


some tolerable level as the background clutter changes while accepting the resulting 
change in signal detection probability. 7 

As necessary as CFAR algorithms are in todays multi-function radars, a price must 
be paid for the control of the false alarm rate. CFAR algorithms reduce the radars 
probability of detection by increasing signal-to-noise ratio (SNR) requirements, and can 
also severely degrade a systems range resolution capabilities. One of the more serious 
drawbacks of CFAR processing is its inherent susceptibility to various types of Elec- 


tronic Countermeasures (ECM). 


B. PRINCIPAL CONTRIBUTIONS 

The principal contributions of this thesis lie in the area of CFAR signal processing. 
The first contribution is a presentation of the weaknesses and vulnerabilities of the most 
popular CFAR architectures on a common ground (including jamming). The probabil- 
itv of false alarm and probability of detection performance for each CFAR 1s compared 
using both an envelope approximation detector and a square law detector. In this 
manner, the implications of using a less complex detector can be evaluated. The com- 
prehensive knowledge derived from this investigation 1s then applied as a stepping stone 
toward the design of an original CFAR system suggested by the author. This new 
CrAR architecture called excision greatest-of (EXGO) is shown to be robust in all op- 


erating scenarios including ECM. 


C. THESIS OUTLINE 

The intended breadth of this thesis 1s from a clear explanation of the most basic 
CFAR function to the creation of an original and complex CFAR detector. Chapter I] 
begins with a comprehensive glossary of many CFAR related terms used in the subse- 
quent chapters. A description of a generic CFAR device is then presented. This de- 
scription covers the various subfunctions found in almost all CFAR systems. The 
second chapter concludes with a discussion of the statistical techniques used to describe 
the capabilities and performance of CFAR systems. Chapter III introduces the inherent . 
vulnerability of CFAR algorithms. By covering this material early in the thesis, the 
reader can better appreciate the complex tradeoffs required in CFAR design. Chapter 
IV comprehensively covers the eight popular CFAR detectors. The capabilities of each 
detector will be carefully described. Stressed here are the adaptive threshold approach 
to managing the false alarm problem. Other approaches such as clutter mapping and 
non-parametric detectors are also introduced in less detail. The performance of these 


eight CFAR algorithms are compared using both square law and an envelope approxi- 


mation detectors. Chapter V then compares the major CFAR families with the intent 
of finding the best CFAR algorithm for some particular application. The systems are 
compared under three operating conditions. These conditions are homogeneous noise, 
clutter edges and multiple target environments. Chapter VI presents new results relating 
to the creation of an original and robust CFAR system. In this chapter the envelope 
approximation detector system is studied to determine the best coefficients for the ap- 
proximation for the new CFAR system developed. In Chapter VII the new CFAR sys- 
tem operation and performance is discussed. Also, the new detector 1s compared to 


other CFAR architectures showing relative capabilities. 


Il. GENERAL CFAR CONCEPTS 


A. INTRODUCTION 

This chapter is designed so that a reader can easily discern the important elements 
and theory of CFAR algorithms. A comprehensive set of applicable CFAR related ter- 
minology is included. These definitions will serve as a basis of understanding for the 
entire thesis. The chapter will also introduce general approaches that radar engineers 
use when deciding on which CFAR processor best suits their particular need. A simple, 
generic Adaptive Threshold CFAR processor will then be covered. All the important 
elements that adaptive threshold CFAR algorithms employ will be described in detail. 
A clear understanding of these subelements will be critical in determining a CFAR sys- 
tems operating characteristics and capabilities. This chapter will end with a discussion 
of CFAR system performance criteria and design criteria that will be followed when 


analyzing the many CFAR architectures discussed in this thesis. 


B. DEFINITIONS. 
The following definitions form a basic glossary of detection and data processing 
terms. The knowledge of these terms 1s required for an easy understanding of this thesis. 


This glossary was taken in part from Schlehers work [Ref. 1]. 


1. Adaptive Threshold CFAR 

A processor which provides a constant false alarm rate (CFAR) in a varying 
nonhomogeneous clutter and noise interference environments by adaptively ad- 
justing the detection threshold. The procedure assumes that the general form of 
the interference’s probability distribution is Known except for a small number of 
unknown parameters. The unknown parameters are estimated on a cell-to-cell 
basis by examining the reference cells surrounding the cell under test. The resulting 
estimated interference probability distribution function is then used in each test cell 
to obtain a threshold setting that provides the desired false alarm rate. 


2. Automatic Detection Radar 

In an automatic detection radar, the target reports are formed in the radar’s 
signal processor before interaction with a human operator or further data process- 
ing. The function of the automatic detector is to process the high bandwidth raw 
radar return, which is usually contaminated with many forms of interference, and 
to extract the low bandwidth target reports with a minimal number of false reports. 
The automatic detection process places stringent requirements on the design of ra- 
dars incorporating this feature with particular emphasis on the signal processing 
function. 


oF 


Clutter 

Clutter is defined as a conglomeration of unwanted radar echoes. Surface or 
area chitter consists of reflections from distributed surfaces intercepted by the radar 
antennas mainbeam and sidelobes. Discrete clutter is returned from stationary 
Objects such as water towers, buildings and other unwanted fixed targets, and can 
be very large. Volumetric clutter is reflected from weather, Chall, ane -other atnio- 
spheric disturbances. Angel clutter is primarily reflected from birds and insects, and 
can be very disturbing. Clutter is characterized by its mean or median equivalent 
backscattering cross-section, amplitude probability density function and its power 
spectral density or equivalent autocorrelation function. 


CFAR Loss 

CPAR loss is defined as the loss of detection sensitivity caused by the effect 
of the CPAR processor on the signal or the detection threshold. In general, the 
CAR loss is a function of the number of reference cells used to estimate the un- 
known parameters, the design probabilitics of detection and false alarm, the num- 
ber of iitegrated pulses and the probability density function of the interference. 
Figure | (Ref. 2] illustrates the relative decrease in CFAR loss in dB (for a Cell 
Averaging svstem) as the number of reference cells used in the noise estimation 
process increase at various false alarm rates. [Ref. 3] is the impact CFAR loss has 
on SNR requirements. 


CFAR loss (cB) 





Number of reference cells, M 


Figure 1. CFAR Loss 


A nonlinear device followed by a low-pass filter which extracts the envelope 
of the high-frequency carrier. The envelope detector ignores any phase information 
contained in the carnier and hence provides only amplitude video information. 


. False Alarm 


An erroneous radar target detection decision caused by noise or other inter- 
fering signals exceeding the detection threshold. In radar statistical theory, the 
decision domain 1s divided into a decision between the hypothesis of noise or that 
of signal-plus-noise. The decision boundary is formed by the threshold level de- 
termined by the desired false alarm probability. A type I error or false alarm occurs 
when the threshold level is exceeded and the noise-only hypothesis is in effect. 


. Monte Carlo Simulation 


The technique of selecting numbers randomly from one or more probability 
distributions for use in a particular trial or run in a simulation study. The system 
Or process to be studied is represented by a model which defines over time its es- 
sential characteristics. The model may be manipulated in ways impossible or i1m- 
practical to perform on the system being represented. The dynamics of the 
behavior of the system under stucy may be inferred by the operation of the model. 


. Moving Window Detector 


A scanning radar detector which accumulates the last m radar return pulses 
within each range resolution cell. This 1s accomplished computationallv through 
formation of the running sum of # radar pulses bv adding the latest pulse to the 
accumulator while subtracting the pulse which occurred n PRI periods in the past. 


. Target Fluctuation 


Variation in the amplitude of the target signal, caused by the changes in target 
aspect angle, rotation, or vibration of target scattering sources, or changes in radar 
wavelength. 


CFAR APPROACHES 


The overall CFAR process consists of a series of techniques used by automatic de- 


tection radar systems to control the rate at which target declarations are falsely declared. 


It is the statistical nature of the radar background that makes a number of false alarms 


inevitable. These background interferers arise from receiver noise, clutter (land, sea, and 


rain), ECM (chaff and jamming), and interference from neighboring radars. When au- 


tomatic detection is performed in homogeneous and nonhomogeneous clutter and noise 


environments a combination of actions are taken in the radar to Jessen the false alarm 


effect. The hardware subsystems that reduce false alarms include [Ref. 4]: 


e Transmitter: 
Waveform selection and frequency agility. 


e Antenna: 
Control of sidelobe patterns. 


@¢ Receiver: 
Rejection of wideband interference, matched filters, and Sensitivity Time Con- 
mois 1 C) 


e Signal Processor: 
Moving Target Indicator (MTTI), clutter sensors, and doppler filters. 


e Detection and Data Extraction: 
Adaptive Thresholding, Nonparametric and Clutter Map CFAR. 


This thesis will focus on the detection and data extraction systems of the radar. In 
relation to CFAR, these are the most direct methods of controlling the false alarm rate. 
There are three maim approaches to handle the detection and data extraction portion of 
the CFAR problem. They include Adaptive Threshold, Non-parametric Thresholding, 
and Clutter Mapping CFAR. Adaptive Thresholding and Non-parametric detectors are 
based on the assumption that homogeneity exists in range around the cell under test. 
These are therefore spatially significant techniques. The Adaptive Threshold technique 
assumes that the noise background/density is Known except for a few unknown param- 
eters. The neighboring reference cells are then used to estimate the unknown parame- 
ters. A variety of CFAR designs are addressed in this category. They differ in principal 
according to the assumptions which are made regarding the characteristics of the back- 
ground noise and the parameters chosen to satisfy the requirements of 1mplementation. 
These different architectures estimate the mean level of the noise differently stnce the 
target, clutter and noise can take on various temporal and spatial distributions. Non- 
Parametric detectors obtain a regulated false alarm rate by ranking the reference cells. 
Under the hypothesis that all the reference cell samples are independent samples from 
some unknown density function, the test sample has some uniform density function and 
consequently some threshold level which yields CFAR. Clutter Maps, which are 
temporal CFAR systems store average background levels in numerous range-azimuth 
cells. If new updated levels exceed the average background by some specified amount, 
a target is declared in that range-azimuth cell. Each of these three major CFAR ap- 
proaches will be further described later in the thesis with the emphasis on the Adaptive 


Thresholding detectors. 


D. A GENERIC ADAPTIVE THRESHOLD PROCESSOR 

Adaptive Thresholding CFAR detectors share common processor subelements. In 
ligure 2, many of these tmportant parts are shown operating together as part of a ge- 
neric CVAR Adaptive Threshold system. The major subelements include the input en- - 
velope detector (such as linear-law, square-law or envelope approximation) , a sliding 
reference window that covers the leading and lagging sets of reference cells (called 
neighborhood | and 2), the threshold multipher (T), and finally a comparator that 


compares the cell under test (identified as Y) with the system estimate of the noise power 


level. A discussion of these major subelements follow. 
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Figure 2. Generic CFAR Detector 


1. System Input 
A basic problem that CFAR systems overcome is graphically depicted mn Figure 
3. Ilere a sequence of one dimensional samples representing radar pulse return energy 
levels from various ranges or azimuths is shown. The peaks in the data samples indicate 
the possible presence of targets dispersed with background clutter and noise. Informa- 
tion of this type forms a steady stream of data into the CFAR processing unit which 


must determine which peaks are actual targets and which peaks are false alarms. 
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Figure 3. CFAR Data Stream 


2. Input Detectors 
The operation that frequently occurs tn radar signal processing 1s the computa- 
tion of the magnitude of the input complex sample stream. The processing of both the 
m-phase (1) and quadrature (Q) channels generally allows for an exact representation 
of the signal with no loss of sensitivity. Various choices of magnitude detection schemes 
are available to the radar engineer to best manipulate the I and Q information. Tigure 
2 depicts three common envelope detector choices. The square-law envelope detector 


esumates the magnitude as: 
2 
R=P4+0° | 


The Iinear-law envelope detector estimates the magnitude as : 
R=/P+Q°. 2 


Since the digital computation of the square root of the sum of squares of the quadrature 
components is a complex task, many approximations to the envelope detector have 
evolved which are linear combinations of the quadrature components. Filip [Ref. 5] has 
shown there are 13 useful approximations which have been devised to satisfy various 


criteria. The envelope approximation detector estimates the magnitude as: 


R=ax Max{|N,JOl] + 6 x Min{|,1Q1] 3 


where a and b are scaling values. Tables 1 thru 3 describe the output Probability Den- 
sity Functions of noise and signal plus noise for each of the detectors assuming that the 
input noise 1s normally distributed N(0,1).. More information concerning input detectors 


will be covered in Chapter VI. 


Table 1. PDFS FOR A TEST AND REFERENCE CELL FOR A SQUARE LAW 
DETECTOR 
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Table 2. PDFS FOR A TEST AND REFERENCE CELL FOR AN ENVELOPE 
DETECTOR 


Detector Output aan 


Noise Only PDF 


Signal Plus 
Noise PDF ae an AO yt 


2p° 


Ded) = ra CXP 


A = amplitude 
I) = Modified Bessel function of order 0 


B = standard deviation of noise 





Table 3. PDFS FOR A TEST AND REFERENCE CELL FOR AN ENVELOPE 
APPROXIMATION DETECTOR 
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(Table Continued) 
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3. Reference Cells 

In our generic CFAR system, the sampled signals in range are the input to the 
reference cells. These neighboring cells yield an estimate of the true noise level in the 
cell under test (identified as Y). By estimating the nearby noise levels, one can determine 
the amplitude difference between the test cell and its neighboring cells and thereby de- 
termine if a target is present. 

Often the two reference cells that are immediately adjacent to the cell under test 
are considered guard cells. A closer took at the guard cells is shown in Figure 4. These 
cells are often used to ensure that signal energy of the test cell does not spill into the 
adjacent cells and affect the noise power estimate. Guard cell values are therefore 1g- 
nored in many Adaptive Threshold processors. Another question that often arises 1s 
whether or not the test cell itself should be included in the noise estimation. Physically 
speaking, the implementation of the system is often simpler if the cell under test 1s in- 
cluded in the estimate. Unfortunately, the statistical representation and theoretical 
analysis of a CFAR processor is much easier when the test cell 1s excluded. Studies 
however, have shown that the exclusion of the test cell from the reference cell summa- 


tion can be used to accurately predict the performance of a system that physically in- 
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cludes the test cell. A’ slight modification to the system threshold value and proper 


choice of the number of reference cells to average are sufficient for accurate predictions. 
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Figure 4. Reference Guard Cells 


One of the more important considerations in Adaptive Threshold CFAR archi- 
tectures is the appropriate choice of the number of reference cells to use. It is easily 
shown that as the number of cells utilized in the estimate of the mean clutter level in- 
creases, the probability of target detection approaches that of the optimum detector 
where the mean level of the clutter plus noise 1s known a priori |Ref. 6]. A tradeoff 
however must be considered since too many reference cells results in greater signal 
processing time as well as an increase in the probability of entering or crossing a clutter 
cdge region. This is a result of an increased physical size of the reference cell region. 
Also of concern is that the likelihood that an interfering target or a large clutter return 
entering the reference window increases with a larger choice for M (total number of 
cells). Another important consideration is that a high number of reference cell samples 
will result in the inevitable violation of the assumption that the noise samples are iden- 
tically distributed over the entire reference window. A desirable goal 1s to use enough 
reference cells so that the CFAR loss is less than 1 dB, and at the same time not Ict the 
reference cells spatially extend beyond one nautical mile on either side of the cell under 
test [Ref. 7]. 

4. Noise Estimate Calculation 

The calculation of an estimate of the mean noise power (Z) can be accomplished 
by a great varicty of techniques. It is the estimation process itself that differentiates the 
various CFAR processors. The type of processing done here is highly dependent of the 
specific clutter and interference models assumed, particularly in a non-homogencous 
environment. Chapter IV deals specifically with the varieties of clutter and noise esti- 


mation techniques. 
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5. Thresholding 

A desirable CFAR scheme is one whose probability of false alarm is insensitive 
to changes in the noise power within the reference cells. In this case, the optimum de- 
tector sets a fixed threshold to determine the presence of a target under the assumption 
that the total homogeneous noise power mean is known a priori. In reality, the calcu- 
lation of the threshold must not only make an allowance for the specified probability 
of false alarm but also for the varying clutter power in the reference window. In the 
presence of clutter, a fixed threshold results in an enormous number of detections and 
will possibly saturate the data processing capability of an automatic tracking system. 
In fact, any small increase in total noise power can result in a corresponding increase 
of several orders of magnitude in the false alarm probability [Ref. 8]. Figure 5 shows the 
probability of false alarm as a function of the increase in noise power density for a fixed 
threshold. As shown, the false alarm rate increases by a factor 10,000 for only a 3 dB 
increase in the noise power density when a fixed threshold is set. A solution to this 
problem is to use adaptive thresholding techniques (CFAR) that adjust the threshold 
value in the presence of interference to maintain some specified false alarm rate. This 
adjustment 1s accomplished by multiplying the reference cell noise estimate by some 
scaling factor. The value of the scaling factor is carefully chosen by the system engineer 


to manipulate the noise level estimate so that a constant false alarm rate is maintained. 
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Figure 5. Fixed Threshold Loss 
Source: Levanon, N. Radar Principals , pp. 226, John Wiley and Sons, N.Y., 1988 


A graphical display of the thresholding concept is shown in Figures 6a and 6b. 
In 6a, the thresliold value is set at a level above all but the two highest receiver output 
peaks so that any signals that exceed the threshold are assumed to be targets. In this 
scenario we have one missed target and one false alarm. Clearly, a decrease in the 
threshold level will increase the probability of detection but also result in an increase in 
the false alarm rate. In Figure 6b, it is shown that the lowered threshold results m one 
false alarm, yet also vields a 100% detection rate. 
a. Stenal-to-Notse Ratto (SNR) 
SNR is the ratio of average signal power to noise variance or noise power 
[Ref 9]. Although the SNR is not directly inanipulated or controlled by a CFAR sys- 
tem, it plays a very important role in CFAR systein performance. Referring again to 


[igure 6a, it can readily be seen that if the SNR of the input were higher, (implying 
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Figure 6. Adaptive Threshold 


larger signal peaks) then the smaller target return would cross the origmal threshold 
value and result in a higher probability of detection. This fact shows the three way de- 
pendence between SNR, threshold values and detection probabihty. Mathematically, 


a... ; A? 
SNR is the ratio of average signal power (5). to the noise variance or noise power 
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6. Comparator 
This is the final subelement of our Generic CFAR detector. The role of the 
comparator is simplv to determine the relationship between the adaptive threshold level 
and the cell under test (Y) in determining the existence of a target. In expression form 
we have 


Hy, 


BE, 5 
Ag 


The comparator compares the cell under test with the threshold value (TZ) where T is 
a selected threshold multiplier that ensures that design false alarm rates are achieved, 
and Z is the noise estimate derived from the neighboring reference cells. The notation 
H, denotes the presence of a target in the test cell while //, 1s the null hypothesis (noise 


only). 


E. DETECTION STRATEGY AND PERFORMANCE CRITERIA 

Early in the design of a radar system an appropriate performance criteria must be 
created for the systems-unique detection problem. This will serve the radar engineer in 
three ways. First, it sets a focus toward the desired properties of the system, it forces 
the design engineer to set up quantitative measures on which to base the design, and 
finally gives rise to specific detector structures that can be implemented or with which 
other suboptimum schemes can be compared. 

The design of CFAR processor for use in radars that operate in clutter, interference, 
and jamming environments require careful consideration of many factors which will af- 
fect the systems performance. CFAR processor designs are always a compromise be- 
tween hardware complexity, CFAR loss and CFAR performance. Some of the factors 
or limitations used in the development of a detection strategy are : 

e least average CFAR loss 

¢ small SNR ratios for detections at long ranges 
e hardware and software constraints 

® processing speeds 

e information storage capacities 


e ability to contain the false alarm rate in the presence of jamming 
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A simple truth of CFAR is that in high clutter and jamming environments, the funda- 
mental purpose of the detection strategy must be the control of the false alarm rate over 
a wide range of variations; even at the expense of detections. 

Many computation methods and models have been used to obtain performance 
characteristics of CFAR techniques. Usually each method was developed for a specified 
type of fluctuating target. Methods include direct numerical integration, Edgeworth se- 
ries, recursive methods, Monte Carlo techniques, and interpolation based on curve fit- 
tings. When using models to estimate CFAR detector capabilities two main concepts 
must be remembered. First, that the model chosen in designing CFAR detectors will 
significantly affect detector performance, particularly when statistical uncertainty exists. 
Secondly, that it 1s impossible to describe all radar working conditions into a simple 
model which will inevitably lead to system estimation errors, 

For the purpose of this thesis and the examination of the many CFAR algorithms 
available today, three scenarios will be analyzed. The Gaussian noise scenario, the 
clutter edge scenario and the multiple target situation. These scenarios will be used to 
evaluate the performance of the most popular CFAR algorithms on a common ground. 
These conditions represent the three most important environments for the CFAR 
processor and can either occur naturally or be artificially generated by a jammer (ECM). 
The Gaussian noise model describes the situation where the radar is thermally noise 
limited. In such a model there are two interesting cases. The first is when the target 1s 
in the clear but the reference cells have background noise, and the second is when there 
is uniform background noise over the entire reference window (including the test cell). 
In both situations the assumption is made that the cells of the reference window are in- 
dependent and contain the same statistics. The clutter edge model is used to describe 
and study transition areas between regions with very different noise characteristics. 
These transitions occur naturally and can be found throughout the reference window. 
Various distributions including log normal, Weibull, and the K distribution can be used 
to represent clutter edging. Lastly, the multiple target situation occurs occasionally in 
radar signal processing when two or more targets are at a similar range. The consequent 


masking of one target by the other is called suppression. 


F. CFAR STATISTICS - PROBABILITY OF DETECTION AND FALSE ALARM 
The radar target detection process is inherently probabilistic or statistical in nature. 
This is due to actual targets being intermixed with randomly fluctuating noise levels. 


Often it impossible to ascertain if an increase in receiver output is the result of a target 


1) 


appearance or the result of noise activity. It is possible however, to declare probabilities 
for this detection process and establish some quantitative values. As defined earlier, the 
probability of detection (/,) 1s the probability that the signal when present 1s detected. 
The probability that some noise fluctuation wili be mistaken for a target is called the 
probability of false alarm (?,,). These two values form the foundation for CFAR sta- 
tistics and analysis. 

As shown previously in Figure 6(a), the threshold value 1s characterized by a voltage 
V (from Figure 2 V,= Tx Z), which when exceeded results in target declaration. There 
is always a probability that this threshold voltage will be exceeded when no real target 


is present. The probability of false alarm can be found from the equation 


Pram | Pu(t) de 6 
V, 


where p,(v) is the probability density function (PDF) of the noise. The probability of 
detection is given by the similar expression as the PDF 1s that of the signal and noise 


combined 


Pq = | Psnlv) dv. 7 
V, 


t 


Figure 7 details an example of these two PDFs that overlap each other. From this plot 
we can easily see the P,, region where the noise statistics are greater than the threshold 
level. Also it is clear that some actual targets are below the threshold level and are not 
detected giving a P, less than 1.0. The signal plus noise PDF (p,,) depends on the SNR 
as Well as the signal and noise statistics. Thus, the single pulse detection probability can 
be also expressed as a function of the signal to noise ratio. For example, with envelope 
detected input the single pulse probability of detection can be described as : 


V, eens 
ak (SNR)2)}. 8 
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Figure 7. Two Superimposed PDFs 
Source: Levanon, N. Radar Principles , p. 42, John Wiley and Sons, N.Y., 1988 
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IW. CFAR ARCHITECTURE VULNERABILITY 


A. INTRODUCTION 

As necessary as CFAR algorithms are in automatic detection and tracking radar 
systems, a loss is incurred with their use. A price 1s be paid for the elimination of the 
human operator. CFAR algorithms are not the quick fix a casual reader of the literature 
may assume. There are real costs incurred at various levels among the numerous CFAR 
algorithms available today. The causes of CFAR limitations will be discussed in detail 
in this chapter and will form a knowledge base to be referenced when comparing CFAR 
systems. The major problem areas include errant operational environment assumptions 
made by the radar engineer, clutter sources and edges, and the multiple target situation. 
Also, the powerful effect ECM has on CFAR will be covered. Electronic jamming has 
the umque capability to replicate at will, those conditions that naturally plague CFAR 
algorithms. The loss effects incurred by CFAR detectors include decreased SNR re- 
sulting in a loss in system P, In sum, one can claim that the CFAR action may ‘sup- 


press’ many real targets in its quest for false alarm control. 


B. VULNERABILITIES 
1. Operational Assumptions 

The design of the detector for a CFAR system significantly affects CFAR per- 
formance. The general operational assumptions made in most CFAR algorithms 1s that 
noise and clutter energies found in the reference windows fall into one of two categories. 
First 1s that the noise is homogeneous, where the statistical parameters of each cell are 
identical [Ref: 10]. Second, that the interference fields are heterogeneous (fields having 
widely dissimilar elements); yet the functions controlling these parameters are known a 
priori [Ref. 11]. Clearly, in actual operation these conditions are violated; that 1s, the 
interference environment is mismatched yielding interference statistics that differ from 
the assumed model. In this case it can be shown that system performance significantly 
degrades, the false alarm rate is no longer maintained and serious target masking may 
be introduced. Another general assumption is that no other signals except receiver noise 
are present in the neighboring reference cells. Clearly this assumption 1s often violated 
as Well in any dense operating environment. Also many types of clutter, particularly 
clutter with high specular reflectors or only a few dominant scatters have been known 


to cause significant deviation from model values. The presence of these strong signals 





in the reference cells has a serious affect on system performance particularly if the 
interferers are stronger than the desired signal. 
2. Clutter and Edging 

Background reflectors, undesirable as they are from the standpoint of detection 
and tracking, are generally denoted by the term clutter. Clutter, which tends to occur 
in contiguous patches forms the basis for non-homogeneity in the system background. 
Continuously distributed in the form of a rain cloud, or manifested as spikes in individ- 
ual cells. Non-uniform, strong clutter is one of the most severe problems for CFAR al- 
gorithms [Ref. 12]. Clutter is generally comprised of a continuum of scatterers, from 
discrete quasi-specular to distributed and diffuse [Ref. 13]. Discrete sources are water 
towers, buildings, sea waves and small hills. Distributed sources are sea echo and rain. 
Table 4 gives four popular density functions for characterizing clutter. Figure 8 displays 
these clutter areas and highlights the differing statistical distribution each contains. The 
figure displays the space -time characterization of clutter data which can be scen as ci- 
ther a succession of spatial snapshots or a bundle of temporal sequences. It 1s worth 
noting that temporal distribution and correlation are relevant for the temporal CFAR 
approach (Clutter Maps) and that the corresponding spatial characteristics apply to the 
spatial thresholding approaches (Mean Level Detector CFAR). 

The clutter edging concept is the effect where within a very small range interval 
clutter levels vary drastically. As shown by Nathanson [Ref. 14] rain, a common clutter 
source, can change intensity from approximately 4 cm/‘hour to over 16 cm/hour and 
back to 4 cm/hour within a one mile range interval. This rate of change of results in a 


clutter return power as high as 60 dB/mile making for a very large clutter edge. 
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Table 4. FOUR POPULAR CLUTTER DISTRIBUTIONS 


Distribution Expression 


Ravleigh y 2 
p(A) =—exp— Az>0 
ae 2A6 
A = amplitude 


As = average RCS 


Weibull 


A = signal amplitude 
B= scaling parameter 
| C= slope parameter 


Log Normal 


] = (ina 


p{A) = ex 


pAJ2n 2p° 
a = mean 

B = standard deviation 

A = amplitude 


K-Distribution 


A = signal amplitude 
Ke_; = amodified Bessel function 
B= scaling parameter 


C = slope parameter 
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Figure 8. Clutter Types 


Two clear effects of clutter and edging on CFAR systems are apparent. First, 
if the cell under test is in the clear while the reference cells are nnmersed tn clutter, a 
masking effect results. That is, the adaptive threshold level 1s increased unnecessarily 
and therefore the P?, along with the P, are reduced significantly. ‘This will occur even 
though there may be a high SNR in the cell of interest. The second case occurs when 
thie reference cells are in the clear but the test cell 1s immersed in the clutter. [In this 
condition, the 7, will mcrease. [Ref. 15] 

In terms of Naval radar operation, one niust recognize that sea clutter 1s inher- 
ently different from land clutter in two basic ways. First, the temporal variations of the 
clutter tend to be larger in magnitude, and second the spatial variations tend to be 
smnaller. A smooth sea forward-scatters incident energy so that little is reflected back to 
the radar, Ilowever, the likelihood of having a glassy sea at any given time 1s quite 
small, so backscattered energy is normally received. Levels of backscattered energy are 
directly linked to wave action and is a function of surface roughness. Radar backscatter 
can occur from the sides of waves as well as the small facets superimposed on the waves. 


These are a function of wind. 
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Numerous models detailing rain, land and sea clutter levels have been devised. 
These models can serve the radar engineer in determining optimal CFAR system char- 
acteristics based on anticipated geographical environmental conditions. [Ref. 9] 
3. Multiple Target Situations 
Closely separated targets are probable in military operating environments. A 
dense target situation occurs whenever two targets come close in range and azimuth even 
when clearly separated by elevation. In this common scenario, both the target returns 
will be contained in one reference window and possibly even in a single reference cell 
[Ref. 16}. If two targets are in the same cell, they are unresolved and act like one target. 
This situation may lead to two undesirable effects. First, if both returns are co-located 
in the same reference cell, some CFAR algorithms would reject that cell as a large clutter 
return and effectively reject both targets. A second and more common effect of inter- 
fering targets in the reference cells is the erroneous behavior of the adaptive threshold 
level. This is due to the interfering target increasing the systems adaptive threshold since 
the interferences are assumed to a legitimate noise samples. In many CFAR applica- 
tions the presence of a strong return among the reference cells can cause a drastic re- 
GuCTION IM system, | ell) |: 
4. Jamming 
a. Basics 

The basic purpose of Electronic Countermeasures is to introduce signals 
into an enemies electronic systems which degrade the performance of the system so that 
it is unable to carry out its intended mission. Certain forms of ECM or jamming tech- 
niques are uniquely devastating to some CFAR signal processing algorithms and clearly 
abuse their weaknesses to unanticipated noise and multiple target situations. 

There are two fundamental ways to introduce jamming energy into a radar 
system. First the receiver noise level can be raised through the injection of external noise 
through the radars antenna. This jamming can be entered into either the radar’s antenna 
mainlobe or sidelobe. This jamming effects Adaptive Threshold CFAR systems by in- 
creasing the voltage threshold. Noise jamming has the effect of obscuring the radar 
target by effectively immersing it in noise. The second and more complex jamming 
technique forces spurious signals into the radars mainlobe or sidelobes to confuse or 
deceive the system. This has the effect of introducing false targets into the reference cells 
simulating one or more interferers. In both cases, the previously discussed vulnerabili- 


ties to clutter and multiple targets are the clear aim of the jammers. [Ref. 18] 
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b. Noise Jammers 

Basic noise jammers can be broken down into different categories such as 
spot, obscuration or broadband jammers. Spot jamming occurs when only a small 
bandwidth is covered by the jammer, whereas the broadband jammer dilutes the power 
density over a large bandwidth. This enables a greater portion of the electromagnetic 
spectrum to be jammed at the expense of effective radiated power (ERP). 

The effect of noise jamming on radars that use CFAR adaptive thresholding 
is the reduction of the detection probability while the system maintains the preset false 
alarm rate. This effect on the P, 1s the result of a grossly inflated adaptive threshold level 
influenced by the jammer noise energies. This degrades the radar system performance. 

c. False Target Jamming 

A comprehensive coverage of various False Target Generating (FTG) tech- 
niques will not be covered here. The many differences among the FTG systems makes 
it difficult to describe them beyond some simple generalizations. In a generic sense, the 
use of FTGs creates transitory false targets that quickly appear and disappear at seem- 
ingly random ranges and angles of arrival. With the advent of sophisticated smart 
jamming systems (those with Electronic Support Measurement (ESM) systems inte- 
grated with the ECM units) the FTGs are capable of repeating ideal waveforms at exact 
radar pulse repetition intervals (PRI). In this scenario, not onlv are the false targets 
more realistic but the jamming may be more effective to an unprepared radar since it 
uses the non-coherent integration gain of the radar to increase the jamming effectiveness 
fer. 19). 

5. Discussion 

Improvement in target detection brought about through clutter suppression and 
Electronic Counter-Countermeasures (ECCM) can be effected through technical ad- 
vances or the removal of errant assumptions by the radar engineer. As Figure 8 shows, 
the probabilistic models of clutter amplitudes change with environmental conditions. 
The characterization and understanding of radar clutter and its effect on performance 1s 
absolutely essential if the radar designer is able to accurately predict expected svstem 
performance. Therefore, in order to make proper model assumptions it will be necessary 
to identify differing clutter types and to be able to describe then properly such as type, 
size, borders, power and spectral features [Ref. 16]. The attempt must be to understand 
the operating environment of a particular system instead of simply trying to suppress 
undesirable energy returns. Unfortunately, the evolution of CFAR algorithms has not 


taken this approach. Rather new algorithms take advantage of more efficient technol- 


14s 


taken this approach. Rather new algorithms take advantage of more efficient technol- 
ogies and techniques to incrementally improve CFAR performance. These improve- 
ments, which are sometimes quite effective have the cost of additional system 


complexity, processing time, and cost. 
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IV. CFAR ARCHITECTURES 


A. INTRODUCTION 

As indicated earlier, each CFAR system has its own characteristic method for esti- 
mating the noise as well as 1ts own method for determining the adaptive threshold level 
based on that estimate. This chapter describes eight of the most popular CFAR archi- 
tectures. Information on how the systems operate, how they determine their threshold 
levels, their performance plots and a discussion of their inherent weaknesses and 
strengths will be detailed. The adaptive threshold systems will be broken down into two 
major categories: Mean Level Detectors and Ranked Order Detectors. Also included in 
this chapter will be a discussion of Non-Parametric and Clutter Mapping CFAR tech- 


niques. These are non-adaptive threshold systems that have merit and deserve attention. 


B. MEAN LEVEL CFAR PROCESSORS 

In this section three mean level algorithms are discussed. They are the Cell- 
Averaged (CA), the Greatest-Of (GO), and Smallest-Of (SO) CFAR systems The deri- 
vation of the probability of false alarm and probability of detection using square-law 
detection are given for each processor. Probability of false alarm and detection plots 
however, have created via Monte Carlo simulations for these CFAR architectures. Both 
square-law and envelope approximation (with a= 1 and b=1) detector results are com- 
puted for comparison. 

1. Cell Averaging CFAR 

a. Background 
The CA CFAR method as first introduced by Finn and Johnson [Ref. 20] 

in 1968 is the most basic adaptive threshold CFAR algorithm. This system can be 
viewed as the first step in a long evolutionary chain of CFAR systems. The CA CFAR 
method uses the maximum likelihood estimate of the noise power to set the adaptive 
threshold under the assumption that the output of the reference cells are statistically 
independent and identically distributed (IID) random variables. When the operating 
conditions of the radar meet this criteria, the CA CFAR detector is optimal in the sense 
that the P, approaches that of the ideal Neyman-Pearson detector as the number of 


reference cells becomes large [Ref. 21]. 
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b. System Description 

The schematic diagram shown in Figure 9 outlines the CA system, which is 
similar to the generic CFAR depiction previously shown in Figure 2. This description 
shows the CA summations of the left (leading) and right (lagging) reference cells. This 
summing and normalizing action is what makes this detector a member of the mean level 
family of estimators. The summation and averaging is this systems unique way to 
measure the mean noise level. In operation, the returns from a given pulse are detected 
and a sample 1s taken from each range resolution cell. The cell under test is the central 
cell. In the CA CFAR, the inputs of the M number of cells are summed resulting in an 
estimate of the background noise. The adaptive threshold level is obtained by multi- 
plying the summed value by a scaling factor (threshold multiplier) depicted as a. This 
value is then normalized by M yielding the overall adaptive threshold level. The mag- 
nitude of the test cell will then be compared to this adaptive threshold in order to de- 
termine the presence a target. 

c. Statistics and Performance 

As CFAR systems evolve technically, their characteristic statistical repres- 
entations often increase in complexity and length. The statistical representation of the 
CA CFAR Is one of the more simpler descriptions and therefore most easily understood. 
With this in mind a comprehensive examination of CA CFAR statistics will aid in the 
understanding of how other systems are statistically represented. The following 
equations are taken from Levanon’s CA CFAR discussion [Ref. 2]. 

Beginning with the assumption that the M samples from the reference cells 
are independent and Gaussian, the envelope r of Gaussian noise will have the PDF given 


aS 


2 


ip 
el A eae 9 


where f is the noise RMS value and A 1s the target amplitude. By normalizing the en- 
velope with respect to f and accounting for the square law detector by making the 


transformation 


z= 10 
the k,, sample of the normalized detected noise thus has the PDF of 
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Figure 9.) CA CFAR Schematic 


P(z,) = exp' 74), 1] 


Using v for the output of the sumunation of leading and lagging cells 


M 
y= > ze 12 
k=] 


the PDI of y is given as 


t)) 
jai ie 13 
= eye 4 
The threshold is then set at 
7. , 
=o 14 
a at 


where T is the threshold multiplier that determines the probability of false alarm. 


To find system P,, (the probability that the magnitude of the test cell with 
a target will surpass I”,) the Rayleigh PDF of the target amplitude is given as 
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A =A? 
p\A)= ve EXP 24? Is 
0 


where A, is the most probable amplitude relative to the average signal power. The PDF 
of the signal plus noise is given as 
] V4 


p(z)=————_ exp? 16 
1+ (S$) 


M+ (9) f° 





A? 
where —_ is the average SNR. Using (16) and the threshold fixed at V, the probability 


of detection will be obtained as 


ona) = 
PLOY aA, -| p(z) dz = exp! TFSNR ?, 17 
V, 


T1owever, in the CFAR system the threshold 1s a function of the random variable y given 
in (14). Thus (17) is only a conditional probability of detection and the overall P, will 


be obtained by averaging (17) over all y, that 1s 


I PLYAES9) -| PASNR|V,= T= ) py) ay. 18 
y=0 : 


Using (13), (17), and (18) we get 


ae (-T+) oe 
P (SNR) = | exp! 1t+SNR ) (Mi)! exp? dy. 19 
pall 


Finally, this known integral can be represented as 


, o if ~M 
PA(SNR, T, M=(1 4+ MU + SNR) ‘aa 20 


The P,, can be easily obtained from (20) by setting the average SNR to zero. That is 


f ere 
Pra (T, M) = (1 +7) um 21 


a2 


The respective Monte Carlo CA CFAR P,, and P, plots are shown in Fig- 
ures 10 and II for both the square law and envelope approximation detectors. As 
shown, the P,, for the envelope approximation CA CFAR systems has a false alarm rate 
10-? at a threshold multipher value of approximately 2.6, and a 10-4 rate at approxi- 
mately 3.9. The square law detector has much higher threshold multiphers. To main- 
tain the false alarm rates of 10-? and 10-4 the threshold multiphers of approximately 5.0 
and 10.8 are required. Using these threshold multiplier values the P, curves fora M = 
32 cell system are generated. As expected, the superior false alarm rate threshold values 
require significantly higher SNRs to maintain constant detection rates. In the envelope 
approximation system for example, with a false alarm rate of 10-? a P, of 0.6 requires 
approximately 7 dB SNR whereas a 10-4 system requires almost 11 dB. The square law 
systein shows slightly better performance requiring approximately 0.5 to 1.0 dB less SNR 


to achieve comparable detection rates. 


d. Strengths and Limitations 

When noise or clutter 1s stationary in the reference cells, CA CFAR detec- 
tors maintain effective CFAR action. Under these conditions CA CFAR 1s the preferred 
detector im that it optimizes the tradeoffs between the P,, and the P,. When these con- 
ditions are not met the performance of the system decreases as the input interference 
departs from the assumed Rayleigh distribution. 

The two basic limitations of the system stem from interference found inside 
the test cell or among the reference cells. As previously mentioned, when the test cell 
Is immersed in strong clutter regions, and the reference cells are in the clear, a natural 
reduction in the threshold will occur resulting in an increase in the P,. This common 
situation lead to the creation of the Greatest-Of (GO) CFAR algorithm, introduced in 
1972 by Hansen [Ref. 22]. For the second case, when the test cell is in the clear and the 
reference cells have interferers, the threshold will be unnecessarily increased reducing the 
system P,. This led to the introduction of the Smallest-Of (SO) CFAR algonthm intro- 
aueed by Trunk |Ref. 23}. 

2. Greatest-Of CFAR 
a. Background 

As stated, the GO CFAR was developed to overcome a decreased adaptive 
threshold level in response to clutter regions. The GO CFAR attempts to correct this 
weakness by independently measuring background noise levels from the leading and 


lagging reference cells and then selecting the larger of these two values for use in the 


33 


PROBABILITY OF FALSE ALARM CA CFAR 
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Figure 10. CA CFAR Probability of False Alarm 
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PROBABILITY OF DETECTION FOR CA CFAR 
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Figure 11. CA CFAR Probability of Detection 
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Adaptive Threshold process. This process results ina selection of the reference window 
that covers the clutter region thus resulting in a closer estimate of the true clutter envi- 
ronment. 

Of note, GO CFAR has also been termed Max Mean Level Detector 
(MX-MLD) CFAR by Ritcey [Ref. 24]. This naming denotes both the type of local es- 
timator (MLD) as well as the combining operation (Max vs. GO). Both terms are ac- 
cepted. 

b. System Description 

A conventional GO processor is depicted in Figure 12. As shown, the 
square law or envelope approximation detector output is fed into the reference windows. 
The leading and lagging window summations are denoted by Y, and Y, respectfully. 
Both reference cell neighborhoods contain M/2 total cells. The detector threshold volt- 
age level (V’,) is obtained by selecting the greater of Y, or Y, (normalized by M/2) and 
then multiplied by the threshold multiplier T. A target is declared when the cell under 
test exceeds V,. In GO CFAR, the background noise is assumed to be Gaussian, the 
target in the test cell and any interfering targets are assumed to be fluctuating inde- 
pendently, each according to Rayleigh PDFs [Ref 17]. 

c. Statistics and Performance 

The following Square Law statistical description of the GO CFAR has been 
analyzed by Hansen [Ref. 25]. The detection performance of this process is derived as 
follows : 

For a SNR at the input to the Square Law detector, the normalized pdf of 


the signal-plus-noise is 


- | —Xx 
As) =k (1+ SNR)? P+ SNR Jeet. 22 


The noise level estimates Y, and Y, are IID with pdfs 


M exp” 
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and the cumulative distribution function (cdf) of Z is then given by 


F(z) = F(z) Fy (2), 24 


where F(z) is the cdf corresponding to p,(v). Thus 
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Vieuwre 12. GO CFAR Schematic 
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The probability of detection is then 


P= | | pAz) x DASIEXAZ,. 26 
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By direct evaluation of the inner mtegral and termwise integration the followmg numer- 


ical result is obtained. 
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When the SNR 1s set to zero, this expression yields the system P,, as 
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The P,, and P, curves for the GO CFAR is shown in Figures 13 and 14 for 
the square-law and envelope approximation detectors. As with the CA CFAR, these 
plots were generated via Monte Carlo simulation. As Figure 13 shows, the P,, values 
of 10-? and 10-4 yield threshold multipliers of approximately 2.4 and 3.65 for the envel- 
Ope approximation detector and 4.5 and 9.8 for the square law system. These multipliers 
in turn generate the P, plots shown in Figure 14. As always, the superior false alarm 
rate systems require higher SNR to achieve comparable detection rates. For example, 
a 0.6 detection probability requires approximately 7.5 dB SNR at 10? ( P, ), and ap- 
proximately 11.0 dB at 10-* for the envelope approximation system. The square law 
results again show an approximate 0.5 to 1.0 dB improvement over the envelope ap- 
proximation detectors. In the following chapter these values will be compared with 
other CFAR architectures. 

d. Strengths and Limitations 

A key advantage of using the MAX/GO family of detectors 1s that near the 
edge of clutter regions the ‘greatest of reference cells capture the desired clutter samples 
and maintains the false alarm rate. Also when a GO detector operates in a benign en- 
vironment suitable for (CA CFAR), only a small CFAR loss of 0.1 to 0.3 dB is noted. 
This loss is due to the reduction of the total number of reference cells (by half) available 
for noise estimation. 

Unfortunately, GO CFAR maintains the false alarm rate in clutter regions 


at the expense of the multiple target scenario. In the analysis by Weiss [Ref. 15 ]aGO 
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PROBABILITY OF FALSE ALARM GO CFAR 
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Figure 13. GO CFAR Probability of False Alarm 
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Figure 14. GO CFAR Probability of Detection 
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system in the presence of an interfering target shows that target detection is nearly in- 
hibited. This problem can be grasped intuitively since when the interfering target is 
large, the half portion of the reference cells containing this extraneous target 1s almost 
always selected in the GO process, hence the threshold is further increased making for 
greater detectability loss. 
3. Smallest Of CFAR 
a. Background 
The Smallest Of (SO) CFAR system introduced by Trunk [Ref. 23] is de- 
signed to handle closely separated target situations. Generally, if the targets are close 
together, the detection from both targets are inerged and a single target 1s reported. The 
problem of resolving merged targets is not only a function of the target separation but 
also a function of the signal strength. Usually it is assumed that if targets are large 
enough to be detected, then they can be resolved if thev are separated by at least one 
pulse width (PW) or equivalently, lie within different range cells [Ref. 23]. 
b. System Description 
The operation of the SO system Is exactly the same as the GO except that 
the smaller of the normalized leading and lagging reference window summations are used 
as the noise power estimate. Figure 15 shows the SO processor architecture with the 
sole change being the selection logic. 
c. Statistics and Performance 
In the SO CFAR scheme, the noise estimate uses the smaller of the sums 


Peenea y,. That is 


Z=min(Y, ¥) 29 
Where Y, and Y, are defined as 
n N 
ay Ga = oe 30 
I i=n+1 
Gandhi [Ref. 26] completes this analysis stating that the pdf of Z 1s given by 
f2(2) =fi(Z) tA) — V@A@ +A@A)). 31 


Where f, and F, are the pdf and cdf of the random variable Y,. This yields the false alarm 
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Figure 15. SO CFAR Schematic 
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Where ’,,(7) and ‘’,,(7) are the moment generating functions (mgf) of the random 


variables Y, and Y,, T is the threshold multipher, and yw is the background total noise 


power. The detection probability is then obtained by replacing T with T(1+SNR) 


yielding : 
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The SO system performance plots created via Monte Carlo simulation are 
shown in Figures 16 and 17. Figure 16 shows the probability of false alarm using the 
square law detector. As shown, the envelope approximation false alarm rate of 10°? re- 
sults in a threshold multiplier of approximately 2.75, whereas a 10-4 false alarm rate re- 
sults in a multipher of approximately 4.4. The square law system results in multiplier 
values of approximately 2.85 and 13.8 for the 10-? and 10- false alarm rates. Figure 16 
displays the P, curves at the false alarm rates of 10-7, and 10~*. As always, a higher SNR 
is required to maintain the superior false alarm rates. In Figure 17 it is shown that a 
0.6 detection probability (envelope approximation) requires approximately 6 dB SNR 
ata 10-? P,, , approximately 11.5 dB at 10%. The square law results again show a 0.5 to 
1.0 dB improvement over the envelope approximation system. These values will be 
compared against the other MLD systems in the following chapter. 

d. Strengths and Limitations 

The strength of the SO technique is its excellent performance in resolving 
closely spaced targets. The detector performance degrades significantly however if in- 
terfering targets are located in both the leading and lagging windows simultaneously. 
This clearly results in at least one of the interferers influencing the voltage threshold 
value and therefore possibly masking the primary target. Furthermore, the SO processor 
fails to maintain a CFAR at clutter edges. Gandhi [Ref. 26], has shown that a 15 dB 
clutter edge leads to an increase greater than five orders of magnitude in the false alarm 
rate at N= 24 and design P,, of 10-®. Finally, even in a relatively benign environment the 
SO architecture results in an excessive number of false alarms since the SO selection 
yields a very low Adaptive Threshold level. This level is generally lower than many 
clutter spikes, all of which yield false alarms. 

More so than the other Mean Level Detectors, the SO processor 1s highlv 
dependent on the number of reference cells chosen. For a small N, the CFAR loss is 


quite large but decreases considerably for increased N. For example Weiss [Ref. 15] has 
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Figure 16. SO CFAR Probability of False Alarm 
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Figure 17. SO CFAR Probability of Detection 
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shown that the detectability loss of the SO CFAR scheme is approximately 11 dB for the 
NN = 4 case but Only 0-7 tara a 


C. RANK ORDERED CFAR PROCESSORS 
1. Ordered Statistics CFAR 
a. Background 
The Ordered Statistic (OS) CFAR processor was designed to overcome the 
loss in detection performance suffered by CA CFAR when interfering targets were lo- 
cated among the background cells. Introduced by Rohling in 1983 [Ref. 27], OS CFAR 
provides inherent protection against drastic reductions in performance in the presence 
of interfering targets. The OS technique rank orders the background voltages encount- 
ered in the neighborhood areas according to their magnitude and then selects a certain 
predetermined address from this sequence. This address value can be the median, the 
minimum, the maximum or any other value. OS techniques have been proven to work 
satisfactorily in both multiple target and non-uniform clutter areas, although they pres- 
ent a small increment in detection loss. OS CFAR methods overcome many difficulties 
which arise in various situations of multiple targets in clutter, but many detection 
problems in special clutter regions remain to be solved. 
b. System Description 
The schematic of the OS-CFAR system 1s shown 1n Figure 18. In this sys- 
tem the values of the reference cells are first sorted by magnitude. The ordered sequence 


thus achieved is represented by the indices in the parenthesis. 
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Figure 18. OS CFAR Schematic 


The main premise of OS CFAR 1s to use some rank selection from the ordered sequence 
for use as the estimate for the average noise power in the entire reference window. The 
variable A identifies the rank or the address of the cell whose input is selected for further 
processing. ‘Phe threshold level V, is then obtained by multiplying the mput from the 


Aih ranked eell by the scaling factor T so that 


re | 
an 


V,= 77, 


is a random variable and its PDF is a function of the PDI of Z,. The use of statistics 
in OS CFAR processing does not define a single CFAR method but rather a serics of 
several CFAR methods. For any given choice of Z,, a distinct CFAR processor 1s cs- 
tablished. Rohling poimted out that the choice of the representative cell K will effect the 
performance of the OS CFAR without interfering targets. For example, in a 16 cell 
reference window with a single target, the detection rates for a K =10 system are 


markedly different then for a K=14. Generally, with no interferers, performance in- 


creases as K is increased. However, with two interferers, K = 14 has an additional | 
dB loss as compared to the K = 10 case. The poor performance in the presence of 


interferers stems from the fact that we reach the point where M minus the number of 
interferers (J) equals K, implying that the representative cell becomes the highest ordered 


fipect ree reference cell. 
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c. Statistics and Performance 
The statistical representation of rank ordered CFAR becomes more com- 
plex and difficult to grasp as the CFAR architectures expand. The OS CFAR analysis 
is the basis for all other rank ordered systems and 1s shown. Using the previous nota- | 
tion, the PDF of the threshold random variable will be shown, as well as equations for 
P,, and P, The optimum choice of K is the number of reference cells (M) minus the 
anticipated number of interferers (J). The following analysis of OS CFAR was taken 


from Levanons work [Ref. 2] 


K=M-J “(K=209 36 


i — TZx oF 


When Z is a random variable with a pdf p(z) and a distribution function P(z), the Kth 


ranked sample has a PDF 
i - 2 
Px(2z) = K( a [Pe ? il Fe) ae), 38 


where 
a eaweteaey 
De) Sex pe, 39 


and the distribution function is therefore 


Z 
P(z) = | exp’? dz=1-— exp”), 40 
0 


Using 38, 39, and 40 the PDF of the Kth sample is found to be 
A eye 
p,(z) = K( e exp( —z)“—**! [1 — exp? ]®-. Al 
The probability of a noise input in the cell under test crossing the voltage threshold 1s 
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The threshold V, is a function of the random variable z,. Thus the P,, can be derived by 


averaging 42 with V, expressed as a function of z, 


Pa = |e p77”) n(zy) dzp 43 
0 


Or 


Pia = | exp’ *”)p,(z) dz. 44 
0 


Using Equations 4] in 44 the P,, becomes 


T+Al— K\'(K— 1)! 
P, ,~K(4)5 ane kK) ) 45 


K (T+ M)! 


In order to derive the P, we can use the same expression used for the P,, by replacing the 


T with 7” where 


tS mGP ae 
thus 
r= |e pn. (z) dz 47 
0 
or 
1 eek Kees I 
p,=K( It ( JHA — 1) 48 
K (7" + #1)! 
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The P,, for the envelope approximation OS CFAR is shown in Figure 19. 
As shown, a false alarm rate of 10-? 1s found at a threshold multiplier value of approxi- 
mately 2.3 and a 10~ false alarm rate has a multiplier value of approximately 3.6. The 
square law results yield threshold multipliers of approximately 5.0 and 11.3 to ensure the - 
10-? and 10 false alarm rates. Figure 20 displays the resulting P, curves for the envel- 
Ope approximation and square law systems at the false alarm rates of 10-?, and 10-4. 
This figure shows that a 0.6 detection probability requires approximately 8 dB SNR at 
the 10-? rate, and approximately 12 dB is required for the 10~* case. These values will 
be used in the following chapter to compare the different detector tvpes. These Monte 
Carlo solutions are the result of the stimulation where the value of K was set as the 20th 
position of the ordered sequence. 

d. Strengths and Limitations 

In general, the presence of one or more interfering targets among the refer- 
ence cells causes the adaptive thresholds to increase erroneously. In OS CFAR this in- 
crease is relatively small. A reference cell with an input from a strong target will be 
ranked at the top, namely, it will occupy the Mth out of M cells. Thus the interfering 
target effectively reduces the number of reference cells to M-1. In the presence of J 
strong interfering targets, the effective number of reference cells drops to M-J. As stated 
by Levanon [Ref. 28], the detection loss due to the increase in threshold 1s not extensive 
aslongasJ<M-k. 

Unfortunately the OS CFAR system suffers from two main limitations. 
First, small targets are easily missed in the presence of multiple targets. Clearly, the 
chosen K will detect the large targets but will set the V, at a value too high for small or 
distant targets to cross. A second limitation of OS CFAR is its inability to perform at 
clutter edges. Anv sharp clutter edge gives a clear rise in the system false alarm rate with 
respect to the false alarm rate that would be obtained in uniform clutter. 

2. OSGO and OSSO CFAR 
a. Background 

In this section two modified OS CFAR architectures are analyzed. The 
OSGO (Ordered Statistic - Greatest Of) and the OSSO (Ordered Statistic - Smallest Of). 
OSGO CFAR has all the advantages of standard OS CFAR in the nonhomogencous and 
multiple target situations with a negligible additional CFAR loss in the homogeneous 
environment. The OSSO CFAR’s sole advantage is that 1t has the equivalent processing 
speed as an OSGO. Unfortunately, the OSSO can not control the false alarm rate and 


it also behaves poorly in the nonhomogeneous clutter situations [Ref. 29]. 
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Figure 19. OS CFAR Probability of False Alarm 
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Figure 20. OS CFAR Probability of Detection 


2 


b. System Description 
The schematic for the OSGO and OSSO system is 
shown in Figure 21. The OSGO and OSSO algorithm is based upon two assumptions. 
First that the noise estimation reference cells have an exponential PDI’; and second that, 
noise in the cells are independent and homogeneously spread. The OSGO CFAR algo- 
rithm consists of taking the greater value of the two samples (leading reference cell is 
K,, lagging reference cell is K,). Obtained from Order Statistic techniques applied to the 


two neighborhood regions independently. The random variable Z is therefore found as 


Z = max(A, A). 49 


For an OSSO system the algorithm takes the smaller value of the two representative cells 


so that 


Z = min(K, Ky). 50 


From this point the OSGO/SO svstems perform their operations identical to that of the 
standard OS system. 
ce. Statistics and Performance 
Using Rohlings [Ref. 27] expression for the pdf for a kth representative cell 
of a set of M’2 cells, Elias-Fuste [Ref 29] analyzed the OSGO and OSSO functions 


yielding : 
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Figuie 21. OSGO and OSSO CFAR Schematic 
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Figures 22 thru 25 display the performance characteristics of the OSGO and 
OSSO envelope approximation and square law systems. Figures 22 and 23 show the 
false alarm rate versus threshold multipliers. As shown the OSGO envelope approxi- 
mation system obtains a 10~? value at an approximate 2.2 threshold value and a 10- 
value of approximately 3.3. The OSGO Square Law system shows the 10-? and 10 
rates at multiphers of approximately 4.1 and 9.35. The OSSO Envelope Approximation 
systems reach a 10°? and 10~ rates at the higher values of 2.5 and 4.05. The OSSO 
Square Law system shows the 10-? and 10-4 rates at multipliers of approximately 6.1 and 
14.8. As shown in the previous CFAR systems, the Square Law system yields a superior 
detection rate as compared to its Envelope Approximation counterpart. Figures 24 and 
25 show the OSGO and OSSO detection rates. In these Monte Carlo simulations, the 
GO or SO of the 10th cell (out of 16) was chosen between the leading and lagging or- 
dered sets to be used as the representative noise estimate. 

d. Strengths and Limitations 
Both OSGO and OSSO systems have a key advantage in that they both re- 


duce processing time in half. This is due to two specialized sorting processors working 
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Figure 22, OSGO CFAR Probability of False Alarm 
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Figure 23. OSSO CFAR Probability of False Alarm 
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Figure 24. OSGO CFAR Probability of Detection 
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independently on the leading and lagging reference cells. Figures 26 and 27 show the 
capabilitics of these two systems in a test environment that includes interfering targets 
and clutter edging. Taken from Eltas-Fuste [Ref. 29], this test contains 256 reference 
cells, 2 clutter edges of 30 dB extending from the 30th to 190th cells, and three targets 
with SNR values of 19, 54, and 19 dB. The interferers are located in cells 100, 105 and 
MO. An additronal target of 22 dB is located at position 215 outside the clutter cloud. 
As Pieure 26 (OSGO system) clearly shows the adaptive threshold level (dashed line) 
always maintains a value greater than the noise plus clutter level (even at the edges). 
Also, all four targets are detected as they cross the threshold boundary. In Figure 27 
(OSSO system) it is shown that all four targets are also detected but with an unaccept- 
able false alarm rate due to an mability to handle clutter edge effects. 

The OSGO CPAR appears to be a fine substitution for the standard OS 
CAR since it mamtams the key OS system strengths and reduces processing time. The 
OSSQO svstem on the other hand, with its inability to control the false alarm rate, makes 


ita poor CVAR svstem for many applications. 
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Figure 26. OSGO CFAR Performance in a Test Environment 
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Ficure 27. OSSO CFAR Performance in a Test Environment 


3. Censored Mean Level Detector (CMILD) CFAR 
a. Background 
The CMLD system was proposed by Rickard and Dillard in 1977 [Ref. 6] 
and is a generalization of the traditional CA CFAR detector with modifications that 
provide robust performance tn the multiple target environment. This 1s accomplished 
by censoring a sclect number of input samples from the ordered group. Like the MLD 
family, CMLD obtains its noise information cstimate from netghboring resolution cells. 
This combination of OS and CA concepts uses an average of all but the first (or first and 
second) largest noise reference samples (i.c. the largest inputs are censored from the av- 
craging routine). The result of this combination system is an architecture that has robust 
propertics that offers superior performance tn multiple target environments. 
b. System Description 
The schematic of the CMLD system ts displayed tn Figure 28. As always, 
the reference window may contain noise and or return cchoes from an interfering target. 
The primary target return echo ts observed in the cell under test. The output of the 
reference cells |g], i= 1,2,....! are fed into a ranking device which outputs the samples 


in ascending order according to their magnitude to yield the M ordered samples 


61 


Vay S12) So S Uap. SS 


At this point the largest K (f or 2) samples are censored. The remaining M samples are 
combined to form an estimate of the notse level in a procedure identical to that of CA 
CPAR. The estimation ts then multiplied by a constant T (threshold multiplier) to yreld 


the adaptive threshold against which the cell under test will be compared. 
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Figure 28. CMLD CFAR Schematic 


ce. Statistics and Performance 
The analysis for the CMLD detection and false alarm probabilitics was or- 
iginally completed by Barkat [Ref. 30]. The following lists the CMLD systems false 


alarm and detection probability equations. 


SNR 
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where is the inverse of the noise power. The detection probability uses the moment 


generating function (‘’) where 


2 
SNR)(x + 
Wi (y) = USNR HD e 
ul[O(SNR) + x] 
and 
a | 
6(SNR) = ) SAR ° 58 
+ eee 
2 
The probability of false alarm can be determined by setting the SNR to zero yielding 
Pig = ‘V¥(Tp). 2, 


lgeuites: 27 and 30 display the /, and P, plots for the CMLD system. The 
envelope approximation system has a 10 false alarm rate at the threshold multiplier of 
approximately 2.7, and a 10~ value at approximately 4.1. The square-law system shows 
corresponding values of approximately 5.7 and 11.8. In Figure 30, the Monte Carlo 
mm@memecetection curves show SNR versus detection rates at 10-? and 10* P, . To 
achieve a detection probability of 0.6, a SNR of approximately 7.5 dB is required at the 
10-? rate, approximately 11 dB at 10~ for the envelope approximation system. The 
square law system shows slightlv improved performance over the envelope approxi- 
mation system. In the Monte Carlo simulation that generated these curves, two of the 
larger ordered values were censored. This in turn led to 30 reference cells being used in 
the noise estimation process. 

dad. Strengths and Limitations 

The performance of the CMLD system exhibits only small additional losses 
in the homogeneous environment (as compared to CA CFAR) but was shown to be 
quite robust when a single large interferer is presented into the reference window. The 
major limitation to CMLD is that the number of cells used for noise estimation should 
be equal to the total number available (M) minus the actual number of outlving 
interferers (J) 1n the reference window to ensure superior detector performance. It has 
been shown by Barkat [Ref. 30] that the CMLD performance is seriously degraded if the 
exact number of interferers is not censored. This requires a priori Knowledge. As ex- 


pected , the more interfering targets censored results in poorer performance. It is clear 
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Figure 29. CMLD CFAR Probability of False Alarm 
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Figure 30. CMLD CFAR Probability of Detection 
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the if the number if interferers is unknown, the CMLD will not only lose its robustness 
but also its CFAR properties as well. For example, if the ordered list 1s undercensored, 
then the noise estimate will be contaminated and result in a degraded detector. On the 
other hand, if overcensoring occurs, the noise estimate is underestimated and an exces- | 
sive number of false alarms occur. 
4. Trimmed Mean CFAR 
a. Background 

The Trimmed Mean (TM) CFAR scheme is a generalization of the OS 
scheme in which the noise power is estimated by a linear combination of ordered sam- 
ples. In the TM CFAR processor, a symmetric or asymmetric number of cells are 
trimmed or censored from both the upper and lower ends of the ordered list. The 
threshold is then estimated by forming the sum of the remaining cells. In TM CFAR, 
as in CMLD CFAR, the censoring points are preset. Again, this implies that some a 
priori knowledge about the background environment 1s required to sensor efficiently the 
unwanted samples [Ref. 31]. 

b. System Description 

The schematic of TM CFAR is shown in Figure 31. The TM svstem first 
sorts the outputs of all the reference cells by magnitude. Then, it judiciously censors the 
K, lower and K, higher ordered samples in the reference window irrespective of the actual 
background environment. When A, = K,, the system symmetric and when these values 
differ the system follows asymmetric trimming. The noise level estimate for the cell un- 


der test is then set to be the normalized sum of the uncensored samples, that 1s 


Ne 
V “ 60 
LW K Ae) 
j= K, +1 


Gandhi has shown through testing that the choice of K, plays the critical role in a ITM 
system performance. Table 5 [Ref. 26] displays the effect that different values of A, and 
K, have on the adaptive threshold multiplier level (M =24 and P,, = 10-6) for both 
symmetric and asymmetric systems. As shown, when the A, censoring value is too high, 
system performance suffers. This is clearly shown by the high valued threshold multi- 
pliers required to maintain the false alarm rate. As shown, th TM system performance 


is most interesting When asymmetric trimming is employed. As shown, with a fixed K, 
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Table 5.) SYMMETRIC AND ASYMMETRIC TRINMING EFFECTS 
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K, is increased. On the other hand, only minimum degradation occurs when 4K, 1s fixed 


and A, iS increased. [Ref. 26] 
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c. Statistics and Performance 
The following analysis of TM CFAR was taken from Gandhi's work [Ref. 
26]. Beginning with the estimate Z given by 


NZ Rok, 
= os 61 


i=] 


The mgf (moment generating function) of Z is therefore the product of the individual 


megf of the x, ‘s. Therefore the false alarm rate is found to be 





Eee: 
Pa= || ¥(D 62 
i=1 
where 
Y(t, 63 
t a; + i 


where a, = (N- K,-1 + 1)/(N-K, - K,-1 + 1). The detection probability 7 iguoe 
tained by replacing the threshold multiplier T with T/(1+ SNR). 

Figures 32 and 33 display the performance plots of the TM CFAR system. 
The probability of false alarm and detection data were Monte Carlo simulated. As 
shown in Figure 32, a threshold multiplier of approximately 2.6 is required to maintain 
a 10-? false alarm rate and a value of 4.0 1s necessary to maintain a 10-4 rate for an en- 
velope approximation system. The square law results that a threshold multiplier value 
of 5.6 and 12.3 are required to maintain the false alarm rates of 10-? and 10-* The de- 
tection plots shown in Figure 33 show SNR versus detection probability at 10-2, and 
10-° P,,. As shown, a 0.6 probability of detection can be found at a SNR value of ap- 
proximately 7.5, and 11.5 for the two respective false alarm rates. In the square law 
system a 0.6 detection rate can be found at corresponding SNR values of approximately 
7.0 and 11.0 dB. These values are very similar to the CMLD system. They were gen- 
erated by symmetrically trimming the two largest and two smallest ordered reference 
cells. This resulted in 28 remaining cells to be used for system noise estimation. 

d. Strengths and Limitations 
A main limitation of the TM system is that as trimming increases, both the 


scaling factor T and the CFAR loss increase. To compound this problem, a priori in- 
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Figure 32. TM CFAR Probability of False Alarm 
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formation 1s critical in order to properly choose trim rates. Consider the situation where 
the leading half of the reference window contains cells from clutter plus noise and the 
lagging half from a clear background. The noise power estimate will include both clear 
region and clutter plus noise background regions. The corresponding threshold will then 
not be high enough to regulate the false alarm rate if the test cell contains a return from 
the clutter plus noise area. Also in a non-uniform background, a compromise must be 
made in determining the proper trimming parameters. In order for the process to be less 
sensitive to interfering targets, K2 should be set to a value greater than zero, and the 
value of K1 should be small in order to attain good detection performance in the ho- 
mogeneous background. If the concern however 1s to handle clutter edges, K1 should 
be large and K2 should be small. Unfortunately in most cases we are interested in reg- 
ulating the false alarm rate in both clutter edge and multiple target environments. This 
balancing act 1s near impossible to maintain in any dynamic radar operating environ- 


ments making TM performance quite variable. 


D. NON-ADAPTIVE THRESHOLD CFAR TECHNIQUES 
1. Introduction 
Although the focus of this thesis is on the optimization of adaptive threshold 
CFAR architectures, other important systems that perform CFAR functions also exist. 
Discussed here are two such systems. The first technique to be discussed is Clutter 
Mapping (which uses temporal rather than spatial information to control the false alarm 
rate). The second technique 1s a Non-Parametric CFAR that 1s generally insensitive to 
environmental background changes. Although these svstems will not be covered in great 
depth, they are introduced to the reader for completeness. 
2. Clutter Mapping 
a. Background 
The Clutter Mapping (CM) system is a specialized CFAR device which av- 
erages radar returns temporallv over several scans to form an estimate of the mean 
background noise levels. This system 1s quite different from the adaptive schemes that 
use spatially differing inputs. The CM CFAR device compares present returns in each 
cell to a background estimate for that specific cell based on past inputs from that cell 
only. Past inputs refer to previous radar scans. For effective operation, the CM CFAR 
device requires temporal rather than spatial stationarity and is thus ideally suited for 
CFAR problems over land [Ref. 32]. 


71 


Hue to the typically large amount of cells to be processed, the required 
processing specd and high cost of digital memory, early clutter maps were typically of 
the blanking type. This system operated with a simple counter that measured how often 
clutter appeared in a specified cell. When this counter reached a preset level, all returns 
in that cell region Were blanked. Today, reduced costs of digital memory and_ technical 
advancements in signal processing equipment has resulted in a recharged interest in high 
resolution CM CHAR processing. 

b. System Description 

The CM technique sets mndependent threshold levels in each map cell to 

vield a CVAR. As shown in Figure 34, [Ref. 32] this technique operates by dividing the 


radar space into cell units. 


Resolution Cells 


Map Cell - consists of 5 
resolution cells 


Spread Map Cell - 27 calls considered 
during map update 





Figure 34. Clutter Map Range and Azimuth Cells 


l:ach cell represents in range and azimuth one resolution cell. As shown, the five igh 
lighted cells represent one clutter map cell. The spread map cells denote the area con- 
sidered during the map update process. Although only five cells compromise a CM cell, 
due to map spreading a total of 27 cells are considered during map update. 

The CM process consists of a number of steps. First, for each CM cell the 
map spreading logic selects for map update the greatest amplitude of all resolution cells 
within the map cell as well as the additional spreading set that borders the map cells. 
This amplitude is then stored in the proper location of the measurement map. The 
process continues with averaging the content of the measurement map with the current 


stored value in the Clutter Map. The average amplitude in each cell is estimated by a 


fice 


low pass, digital filtering of the input data. Target detection declarations then follow the 
CA CFAR logic form. If the CM cell value exceeds that mean background estimate, a 
target is declared. 
c. Design Issues 

The design of a high resolution CM system is not necessarily feasible for 
any combination of radar parameters. Four key parameters are used in the determi- 
nation of a CM design and include; compressed resolution cell size, radar update period, 
number of pulses noncoherently integrated, and detector law employed. According to 
Farina [Ref. 4], CM CFAR designs are based on : 


¢ The geometry of the map: How the surveillance space is divided into cells (i.e. rec- 
tangular or polar). 


e The map building process: The map changing process may be periodic with in- 
creasing and decreasing counters or bv amplitude averaging. Also, assuming a map 
cell larger than the radar cell, the data obtained in one radar cell may be used alone 
or in combination with data from nearby cells. 


e The approach taken: Different action can be decided on the basis of the content 
of the CM. Such actions include blanking zones, switching to different processing 
channels, resetting detection thresholds in each cell, and tracking clutter points. 

d. Strengths and Limitations 

A primary reason for employing a CM CFAR system 1s for its excellent 
interclutter visibility. Interclutter visibility is the ability to detect and track targets in 
Shadow areas where clutter is normally absent. Another favored capability of a CM 
system 1s that it generally provides better detectability of targets in near tangential flight 
paths over clutter regions [Ref. 7]. Finally, CM systems also may be employed to sense 
locations where clutter echoes are too strong to be suppressed by other signal processing 
systems such as doppler filters. 

The major limitation, or drawback of CM CFAR 1s the assumption that 
clutter statistics are temporally stationary (over five to ten scans of the radar). 
Rainstorms, jamming, and other nearby radars thereby cause excess false alarm rates. 
Also large differences in system capability become apparent even with low velocity point 
clutter. At only two Knots velocity, there is a four orders of magnitude difference be- 
tween the false alarm rate as compared to zero knots. Even with a Jand based radar 
where there should be no apparent velocity between the land clutter and a radar, changes 
in the atmospheric index, multipath and systems instabilities all lead to small positional 


shifts resulting in unacceptable false alarms. 
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3. Non-Parametric CFAR 
a. Background 

Throughout the discussion of the adaptive threshold CFAR techniques, an 
emphasis was placed on optimal detectors requiring an essentially complete statistical 
description of the input signals and noise. Thomas [Ref. 33] mentions three compelling 
reasons which lead to consideration of other ‘nonoptimal’ detectors. First, a complete 
Statistical description of the input is rarely available; second, the statistics of the input 
data set may vary with time or may change from one radar application to another. 
Finally, optimal detectors may be too complex or costly to implement. Adaptive sys- 
tems have developed and evolved to meet the first two conditions and can perform in a 
near Optimal sense in a unknown or changing environment by proper adaptation of de- 
tector structure: however, such detectors tend toward greater complexity. Non- 
parametric or distribution-free detectors exhibit insensitivity to the environment rather 
than adapting to it and often exhibit simplicity in implementation. 

In the Non-Parametric (NP) device it 1s assumed that the statistics of the 
interference are unknown. The rational of the approach is to somehow map the un- 
known PDF onto a known one where a fixed threshold produces CFAR [Ref. 4]. This 
technique enables CFAR performance against very broad classes of noise probability 
density functions [Ref. 22]. 

b. System Description 

A wide assortment of NP processing techniques are available and present a 
the practical problem of choosing the proper technique for a particular need. A com- 
mon and simple NP detector obtains a CFAR by order ranking the test cell among the 
reference cells (Ref. 34]. The smallest ranked value receives a rank of zero and the largest 
a rank of N. Under the hypothesis that the samples are independent with unknown 
PDF, the test cell has equal probability of taking on any of the N ranked values from 
zero to N. The rank detector is then constructed to compare the rank of the test cell 
against a preset threshold rank. If the cell under test rank is greater than that of the 
threshold rank, a target is declared. This simplistic system normally incurs a CFAR loss 
of approximately 2 dB but achieves a fixed false alarm rate for any noise density as long 
as the input sample stream remain independent. 

c. Strengths and Limitations 

The insensitivity to environmental noise density changes is the strong suit 

of any NP detector. Clear costs are paid for this action however. Besides the 2 dB loss, 


correlated samples result in a detector inability to maintain CFAR. Also, a large inter- 
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fering target may lead to target suppression, that is, if a large return is found in the ref- 
erence cells, the test cell will not receive the highest ranking possibly resulting in a rank 
below the threshold rank. A final concern 1s that by maintaining only the rank orders, 
the system loses the actual signal amplitude information which may be used in other 


signal processing applications. 
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V. CFAR ARCHITECTURE COMPARISONS 


A. INTRODUCTION 

A great deal of literature has been written comparing the relative performance of 
various adaptive threshold CFAR systems, each paper dealing with only one to three 
different types. The goal of this chapter is to compare on a larger scale, all eight of the 
popular types of CFAR systems being used today. As always, there are tradeoffs in 
capabilities between systems; that is, even if one detector is superior in one scenario, it 
may be poorer in some other. The system comparisons will be on the basis of capability 
in 1) homogeneous noise, 2) clutter edging, and 3) interfering targets. It should be un- 
derstood that these are but three idealized examples of the multitude of different situ- 
ations which may occur in actual radar operation. 

The comparisons made in this chapter will be restricted to the adaptive threshold 
systems previously discussed. The chapter begins with a comparison of the MLD family 


of detectors followed by a comparison of the rank ordered systems. 


B. MEAN LEVEL DETECTOR COMPARISONS 
The relative capabilities of the MLDs are tested under the three test environments 
of homogeneous noise, clutter edges, and multiple targets. As discussed, the CA CFAR 
is the optimum system in the homogenous noise environment. The modifications to CA 
CFAR (GO and SO) have been proposed to overcome the problems associated with the 
non-homogenous background. The GO system was designed to regulate the false alarm 
rate in the region of clutter transitions and the SO system was designed to resolve two 
closely spaced targets. 
1. Homogeneous Noise 
In the homogeneous environment, the threshold multiplier T can be used to 
judge a systems capability. Table 6 [Ref. 26] details the CA, GO, and SO (Square Law 
systems) threshold multipliers at various false alarm rates and different selections of N. 
The data in this table represents the noise only environment. As shown, the CA system 
has the lowest threshold multiplier in all cases which results in superior P, values. The 
GO svstems values are slightly worse and the SO system values are significantly poorer 
than that of the CA detector. 
Figure 35 details detection probability curves for the MLD family via Monte 


Carlo simulation using an envelope approximation detector. The P, for these curves 
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Table 6. MLD THRESHOLD MULTIPLIERS (SQUARE LAW) 


. N=8 N=8 N=8 | N=16 | N=16 | N=16 
fa (CA) (GO) (SO) (CA) (GO) (SO) 


10-4 2.162 3.600 10.88 0.078 1.360 2.444 
10-6 4.623 LES 36.00 1.371 2.420 oe a 
10-8 9.000 30 117.90 2.162 3.860 202 





was set at 10-* with a total of 32 cells being used. As expected, the CA CFAR is the best 
system in this environment with the GO system showing approximately 0.2 dB additional 
loss and the SO system yielding a 0.7 dB additional loss. 
2. Clutter Edges 

The second comparison area considers detector performance in clutter regions. 
These regions can be caused bv chaff, weather clutter distributed in range, and by patchy 
land clutter. The boundary of this interference (the clutter edge) will move into and out 
of the reference cells as the test cell approaches or leaves the clutter patch. Of prime 
concern 1s the detectors ability to regulate the false alarm rate caused by edging and not 
specifically with detector losses [Ref. 35]. In the troublesome scenario where the clutter 
edge occupies half of the reference cells, Moore [Ref. 35] states that the P,, of the CA 
system increases by a factor of a 1000, whereas the P,, of the GO system only increases 
by a factor of 17.5. This control of the false alarm rate 1s the prime advantage that the 
GO system maintains over CA CFAR. As expected, the SO CFAR system has the 
poorest performance in clutter edge regions. Trunk states [Ref. 23] that the SO 
processor performance worsens by more that 5 orders of magnitude when the clutter to 
noise ratio is greater that 15 dB. Figures 36-38 show the capabilities of the three MLD 
systems confronting a ‘real’ sea environment clutter edge. The adaptive threshold levels 
for each system are shown Versus the actual clutter power levels. Figures 36 and 37 
clearly show that the CA and GO algorithms handle the simulated clutter edge. How- 
ever, the SO system shown in Figure 38 is unable to handle the leading and lagging edges 
of the clutter region leading to unwanted false alarms. 

3. Multiple Targets 

The final scenario to examine is how multiple target situations affect the MLD 

family. The SO CFAR svstem which was designed specifically for the multiple target 


scenario is nearly unaffected by a single interfering target while the suppression 1s serious 
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Figure 35. MLD Family Detection Curves 
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CA CFAR IN CLUTTER EDGING 
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Figure 36. CA CFAR in Clutter Edges 
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GO CFAR IN CLUTTER EDGING 
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Figure 37. GO CFAR in Clutter Edges 
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SO CFAR IN CLUTTER EDGING 
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Fieure 38. SO CFAR in Clutter Edges 
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in the CA CFAR and even worse in the GO. Practically, when the total number of ref- 
erence cells is less than or equal to 16, detection of target pairs with aGO CFAR system 
is almost totally inhibited [Ref. 15]. 

Figures 39 and 40 show the effects of interferers on the envelope approximation 
MLD systems. In Figure 39 a single interferer is introduced into the reference cells. 
As anticipated the SO system has the best performance with a single interferer yielding 
the smallest additional CFAR loss of approximately 1 dB. The CA and GO system 
performances are significantly reduced yielding an additional loss of approximately 3 and 
5 dB respectively. Figure 40 shows the effects of 2 interferers (both in the same reference 
cell neighborhood) on the MLD systems. Again, the SO system maintains its perform- 
ance while the CA and GO systems perform poorly. 

4. Conclusions 

The tradeoffs to be compromised concerning the selection of the appropriate 
type of processor and adequate choice of N are highly dependent on the clutter and in- 
terference models the radar engineer chooses. An optimal and general performance 
CFAR detector can almost never be devised. Therefore it is of great importance to un- 
derstand fully the operating environment of the radar system in order to correctly tailor 
or choose the proper CFAR system that will yield superior results for some particular 


application. 


C. ORDERED STATISTICS VS. MLD COMPARISONS 

As previously discussed, the OS processing scheme was introduced to alleviate the 
problems associated with the MLD family. Of significance is that a properly designed 
OS detector with interfering targets maintains its robustness with only gradual detection 
loss, while suffering only a minor degradation in the homogeneous environment. Like 
the MLD systems though, the OS processor is generally unable to prevent excessive false 
alarm rates at clutter edges unless the clutter appears in one single contiguous patch. 

1. Homogeneous Noise 

In the homogeneous environment the OS CFAR processor performance is 1n- 

ferior to that of both the CA and GO systems. The loss however 1s typically 0.5 dB 
(K= 21, N= 24, P, = 0.5) [Ref. 26] and is quite tolerable. Figure 41 displays the Monte 
Carlo comparison of the CA, GO, SO and OS systems in homogeneous noise. The false 
alarm rate used was 10-4, and the number of reference cells set at M = 32 with K set to 
20 in the OS system. As shown, the OS system 1s only slightly poorer than the CA and 
GO systems but is slightly better than the SO. 
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DETECTION CURVES FOR MLD SYSTEMS WITH A SINGLE INTERFERER 
1 — 
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Figure 39. MLD Probability of Detection with a Single Interferer 
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DETECTION CURVES FOR MLD SYSTEMS WITH TWO INTERFERERS 
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Figure 40. MLD Probability of Detection with Two Interferers 
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MLD and OS PROBABILITY OF DETECTION CURVES 
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Figure 41. Ordered Statistics versus MLD Detectors in Homogeneous Noise 
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2. Clutter Edges 

In considering the clutter edge problem, one would intuitively expect the OS 
CFAR performance to be relatively insensitive to edging. This is true only when the 
clutter returns have slowly varying amplitudes or appear in contiguous patches. The 
ability to handle these edges (such as weather clutter) are detailed in Figure 42. In this 
figure the OS system is able to handle the contiguous patch of sea clutter. Unfortu- 
nately, OS CFAR can not handle random clutter spikes as well. This is caused by clutter 
being found in the cell under test while the representative cell (K) is in the clear. This 
inevitably leads to an increase in the false alarm rate. 

3. Multiple Targets 

In general, the presence of one or more interfering targets among the reference 
cells cause the adaptive threshold to increase erroneously in the MLD family. In OS 
CFAR, this increase is relatively small resulting in a superior performance as compared 
to the MLD processors. As always, the OS processor performance is highly dependent 
on the value chosen for K. An optimum OS system is one where K = M - J, where J 
is the a priori known number of interferers. Figure 43 displays the inherent strengths 
of OS CFAR as compared to the multiple target handling capability of the SO (shown 
to be the best MLD system with interferers) The figure displays probability of detection 
versus SNR for envelope approximation systems. As shown, the OS system 1s favorably 
compared to the SO CFAR system. Both the OS and SO systems can easily handle two 
interferers with little additional CFAR loss. 

To further stress the performance of a properly designed OS system Table 7 
details the detection probability losses and relative CFAR loss of square law OS and CA 
detectors due to J interferers. In this table, both schemes use a false alarm rate of 10-4 
and 64 total reference cells with a SNR of 20 dB for both the primary target and the 
interferers. The value of K is set at 54, allowing for a total of 10 interferers prior to 
system degradation. As expected, with no interferers, the CA system is superior. How- 
ever, once interferers are introduced into the system the loss in probability of detection 
and the relative CFAR loss of the CA system are dramatic as compared to the robust 
OS detector. 

As previously described in Chapter IV, the OSGO and OSSO are direct 
descendents of the OS system. In terms of comparison, OSGO CFAR has all the ad- 
vantages of the standard OS system in nonhomogeneous clutter and in multiple target 
situations with a negligible CFAR loss in the benign environment. Of importance 1s that 


the OSGO system requires only half the processing time that the conventional OS sys- 
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OS CFAR IN CLUTTER EDGING 


e 
OSA 58) 6 Oe Oe Or Oy Oe wise: w) Yeh EPs 6. 8 eee sa, 8) OR 8 eel eee) eo) ere), (Eel ene) ne) we le fe) ele er ce: Dl me a oe in ee) ee a 
e 


Sia). tare el fe) awd fe, wt: Ney Fah ee oe 


SOME Saas ieee eA iO wise) Vet seem MaKe dh, (6LCOn1@ 78 <6" 16.5 (OO SOL ne) ctw ey ielle? tela, ey (0) ee ciate” sisi faiy ede ie i@n ete hwo lee. e@ aes 
e 


SOS CROSS! Oe Me ean ee ONG 8) 8) WNW OS Cet eee 8 Ws. eee 8). fe) “ee, ee ia wi cela | gs ety el @ eel, tec, e) fey le) ao “othe 
Seer eC re on) blaine | is ie) © e, Sime .6 6 ve: @ 6 = 


sile\ us)! ee 18: (eee) ete) s/c ee. eee) a a) a Tes 


Amplitude in dB 


- = = @ = © woe wor ee 
. 


see ce ee Oe Mw et te en AR ee pe et ee te eo Re 


wis) ‘SluelKe, © © (@ife le, e, ©: 6: agi © 10) (8) 6. OO fe ie | (em § neers, fe Up. ve Hes es 008 (ee 6: in 





O 100 200 300 400 
Reference Cells 


Figure 42. OS CFAR in Clutter Edges 
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DETECTION CURVES FOR OS AND SO WITH TWO INTERFERERS 
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Figure 43. OS ys SO with Two Interfering Targets 
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PROBABILITY OF DETECTION CURVES FOR RANK ORDERED FAMILY 
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Figure 44. Rank Ordered System Comparison in Homogeneous Noise 
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Table 7. MULTIPLE TARGET EFFECTS ON CA AND OS CFAR 


a CFAR Loss (dB) CKAR Losi 





tems requires [Ref. 29]. Concerning the OSSO system, its only advantage is that it too 
requires less processing time yet it has much higher loss than the OS system and behaves 
poorly in the non-homogeneous situation. These conditions generally make the OSSO 
system a poor detector choice. Figure 44 shows the P, curves for all the envelope ap- 
proximation rank ordered devices. These Monte Carlo curves were generated at a false 
alarm rate of 10-*. As shown, the OS system with a representative value of K = 20 1s 
only slightly superior to the OSGO system with K = 10. The OSSO system, also with 
K = 10, yields the poorest system performance in this noise only environment. 

From the results obtained so far, a clear conclusion of the OS system perform- 
ance versus the MLD detectors can be drawn. Though the OS CFAR exhibits some 
additional loss of detection in the homogenous noise background, its far superior per- 
formance in multiple target environments makes this a seemingly desirable system. Of 
course, the proper value for K must be chosen to ensure these robust results. This 
generally requires a priori information not generally available, making these systems 


theoretically superior but operationally sub-optimum performers. 


D. CENSORING SCHEME COMPARISON 

The CMLD and TM schemes will be considered together since they are both cen- 
soring schemes that imply some required a priori knowledge to avoid unwanted samples. 
With this knowledge, the TM and CMLD censoring points would result in near equal 
optimal system performance. As previously mentioned, the value of K2 in the TM 
scheme (upper censoring point) plays the crucial role in determining detector perform- 
ance. The non-zero K2 (TM) or K (CMLD) values dictated by robust detector per- 
formance in the multiple target environment conflicts with the requirement to maintain 
the false alarm rate in regions of clutter power transitions. A near zero value for K2 is 


necessary for this case. 
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1. Homogeneous Noise 

In the homogeneous environment, TM and CMLD detectors overall perform- 
ance is better than that of a properly designed OS system and performs nearly as well 
as a MLD system [Ref. 36]. 

Referring back to Figure 44, the CMLD and TM systems show near excellent 
performance in the noise only environment. In these curves the CMLD system censors 
two of the largest ordered cells and the TM censors the two largest and the two smallest. 
As shown in this 10-* case, CMLD and TM outperforms the OS system by approxi- 
mately 0.5 dB and is nearly equal to the optimum CA architecture. 

2. Clutter Edges 

The clutter edge problem for the CMLD and TM system varies directly with the 
censoring points chosen. As previously mentioned, a small censoring point may degrade 
the detection performance since high power clutter samples mask the target in the test 
cell. As always, the inverse problem of over censoring, results in additional CFAR loss. 
Since it 1s recognized that the GO system is the superror MLD system in clutter, it will 
be compared with the two censoring schemes. Himonas [Ref. 21] has shown that in high 
clutter power transition areas TM CFAR with Kl = 0 and k2 = 4 yield almost iden- 
tical performance with that of a GO system when the actual number of clutter cells is 
four. The performance worsens however, as the number of actual clutter cells increases 
or decreases away from the preset values of O and 4. For small clutter power transition 
regions the detection performance of the TM (K1=0 and K2=4) system is actually su- 
perior to that of the GO system by approximately 2 dB. 

Figure 45 shows the envelope approximation TM CFAR system performance 
in the sea clutter edge. Since the TM method sums all but the highest and lowest ranked 
reference cells, its curve is similar to the CA CFAR curve under the same clutter edge 
environment (Figure 36). As shown, the TM censoring scheme handles clutter edges 
quite successfully. 

3. Multiple Targets 

In a multiple target situation, prior knowledge of the number of interfering tar- 
gets will result in superior performance of the TM and CMLD systems. As always, 
system performance decreases rapidly with an improper choice of censoring points. 
Figure 46 clearly proves this point. The envelope approximation CMLD CFAR system 
used in this simulation censors the two highest ranked reference cells. This enables the 
system to handle up to two interferers. As shown, the CMLD system with no interferers 


performs only slightly better than the CMLD system facing two interfering targets. 
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TM CFAR IN CLUTTER EDGING 


, i] 
a ee ee ewe ch hw RM ww ott Oe te 6 Mae Re 8h te 8 ee He 8 ee a et es @ 6 fe 6 6 o 8) Ue) 6 0 0 © © © 8 miss ea = s\ « = «| ss a's) el (allele mealiaemtanee 
P ° 


ry 
Pe Ce ee i OC | 
e 


. 
ove ws eo Ng Se 8 ew 6 8 8 ee wg nl 8 0 8 6 eo 8S fe = ee) es) ie) 0) Sw) ee ee 


re ee) Ce i | 


. 
* = © «© «© «© ote « © oO we he © Fe = F = — s we wo wwe wae 


e O 
cose @ een Cm se it ww 6 eS 8 6 8 8 Ue 8 BR) 4 ee 6 ees eireey 4 @ 0: nm) eee fe a) U=aecae 


"es © ee eee eae 


Amplitude in dB 


aa @ 6 « = =a @ 6 6 6 « 6 4 @ (6156 S18 8) 6 (8&0 ww O80) er 0) (Sie c¢: 8! ©) (a) 19178 10) 0) 6s (O0e a) as Sea ee eS 





0 100 200 300 400 
Reference Cells 


Figure 45. TM CFAR in Clutter Edges 
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However, once the quantity of interferers becomes greater than two, system performance 
suffers greatly. As shown, the CMLD system with four interferers shows a drastic re- 
duction in capability. Figure 47 details the CMLD CFAR loss as a function of the 
number of interferers relative to the optimum Neyman-Pearson detector. In this plot, 
a CMLD system with N = 32, P, of 0.9 and a designed false alarm rate of 10-4 is de- 
tailed. The CFAR loss shown is a result of the increase in false alarm probability caused 
by improper estimates of the actual number if interferers present. For example, if the 
actual number of interferers is four (the x axis) and our system is designed to handle two 
interferers (selecting the curve labeled 2), a CFAR loss of approximately 1.3 dB (read 
off the y axis) occurs. 
4. Conclusions 

The performance of the censoring CFAR systems are robust in all operating 
scenarios as long as the proper censoring points are chosen. When improperly selected, 
system performance degrades intolerably. In actual operating environments, these 
CFAR systems would possibly result in unacceptable detection and false alarm rates due 


to a lack of a priori information. 
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DETECTION CURVES FOR CMLD CFAR WITH TWO AND FOUR INTERFERERS 
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Figure 46. CMLD CFAR with 0, 2, and 4 Interferers 
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VI. ENVELOPE APPROXIMATION RESEARCH 


A. INTRODUCTION 

As discussed earlier in Chapter Two, the input to a CFAR system is often the en- 
velope detected in-phase (1) and quadrature (Q) channels of the baseband signal 
(R=./P+Q’). Since the digital computation of the square root of the sum of squares 
of the quadrature components is complicated and time consuming, various approxi- 
mations to this operation have evolved. One less complex method of I and Q detection 
using absolute values 1s the envelope approximation method where the input is estimated 


as 


R=axmax{|/|,|Q]}+6x min{|/|,|Q]}. 64 


In this calculation, a and b are simple scaling coefficients. 

The purpose of this chapter 1s to examune the performance difference this type of 
detector has on CFAR processors for seven different scaling factors. This examination 
uses a GO CFAR device. Results for the envelope approximation GO CFAR processor 
in terms of probability of false alarm [Ref. 37] and probability of detection [Ref. 38] are 
shown. The a and b multipliers are listed in Table 8 along with the average error and 
mean square error for the seven approximations. In the first five approximations, the 
multiplving coefficients are either one or simply binary fractions. The sixth approxi- 
mation was designed to have zero average error which simultaneously minimizes the 
variance of the error. The last approximation was designed such that the end point error 
equals the absolute value of the peak error in the region 0 < @ < 2z/4 [Ref. 37],[Ref. 


- 


5}. 


B. PROCEDURE AND RESULTS 

Monte Carlo simulations were created to test the seven scalar combinations listed 
in Table 8. The envelope approximation results have been devised for the six cases N 
= 2,4, 8, 16, 32, and 64. The threshold multiphers used in these simulations were taken 
from Paces results to ensure a false alarm rate of 10. Table 9 lists these threshold 
multipliers. Figure 48 shows the P,, curve versus Threshold Multipliers. Shown are the 
seven combinations plotted along with the JP+Q? results. Figures 49-54 detail the 


resulting detection probabilities versus SNR. As detailed in these performance curves, 
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Table 8. SCALING FACTORS 


Case a Scalar b Scalar Average Error |} Mean Square Error 


0 5 
US 


W515 
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the a = 1.0 and b = 1.0 case yields the highest system detection performance and the 
lenvest P,, for a given threshold multiplier. The a = 1.0, b = 0.25 yields the poorest 


performance. 


Table 9. THRESHOLD MULTIPLIERS AT 10-4 PFA 


I ee OAS A CASE 51 CASE 61 CASE? 





C. CONCLUSIONS 


A clear conclusion can be drawn from this study of the envelope approximation 


scaling coefficients. In all cases, a = 1.0 and b = 1.0 yields the best detection per- 
formance. The next chapter uses this information and presents a new CFAR architec- 


ture that overcomes the inherent deficiencies of the architectures examined previously. 
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Figure 50. Envelope Approximation Curves with N = 2 
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PROBABILITY OF DETECTION CURVES - GO CFAR N =4 
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Figure 5t. Envelope Approximation Curves with N = 4 
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PROBABILITY OF DETECTION CURVES - GO CFAR N «8 
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Figure 52. [Envelope Approximation Curves with N = 8 
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PROBABILITY OF DETECTION CURVES - GO CFARN = 16 
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Figure 53. Envelope Approximation Curves with N = 16 
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PROBABILITY OF DETECTION CURVES - GO CFAR N = 32 
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Figure 54. Envelope Approximation Curves with N 
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Vil. EXCISION GREATEST OF (EXGO) CFAR 


A. INTRODUCTION 

As presented throughout this thesis, numerous CFAR architectures have been de- 
veloped to balance the conflicting goals of maintaining a high detection probability while 
enjoying a low false alarm rate. As shown, this difficult task is magnified in the presence 
of interfering targets and ECM. In this chapter a new CFAR device call “Excision 
Greatest Of (EXGO) will be presented. This system was designed to maintaining su- 
perior performance under clutter edge, multiple target and jamming environments. The 
concept of excising large interferers was initially introduced by Goldman and Bar David 


[Ref 39]. for their cell averaging scheme. 


B. SYSTEM DESCRIPTION 

A schematic diagram of the proposed EXGO CFAR detector is shown in Figure 55. 
The EXGO processor uses envelope approximation to detect the inphase (I) and 
quadrature (Q) components of the signal. Scalar values of a = 1 and b = 1 are used 
since they have been shown to most closely match the results of a true envelope detector. 
Two additions to the standard GO MLD are shown in Figure 55. The first addition 1s 
the excision logic, and the second is the extended leading and lagging reference cells 
(shown with hashed lines). 

I. Excision Logic 

The purpose of the excision logic is to compare the relative magnitude of all 

utilized reference cells to an adaptive voltage threshold level (V,,). The excision logic is 
shown in detail in Figure 56. The threshold level is the product of a preset scalar value 


(T,) and the continuously updated running system noise average (R) given by 


pe 65 





where x(j) is the reference window input and k is the total number of inputs into the 


system over time. Thus, the initial threshold level is set as 


V,,=RxT,, 66 
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Figure 55. ENGO Schematic 


If the magnitude of any individual reference cell is greater than V,,, an interfering target 
1s assumed to be present in that cell and therefore ignored from further use in the GO 
process. For example, with 7, = 2, any reference cell that is greater in magnitude than 
twice the svstem noise average will be discarded from further GO processing. If the cell 
is less than twice the running noise level, the reference cell is then processed in a normal 
GO fashion. The selection of a proper threshold multiplier, 7, is important. With too 
low a threshold, proper noise sainples will be excised from the system, thus increasing 
the CFAR loss. Setting 7, too high will cause some large interfering targets to pass the 
excision logic and thus contaminate the noise power estimate and degrade the probabil- 
itv of detection. A binary integration at the output of the excision logic counts the 
number of threshold crossings. 

The EXGO architecture takes special precautions in order to maintain the 
proper false alarm rate in the presence of clutter edges. When at least one entire leading 


or lagging reference cell window is fully contaminated by clutter (recognized by greater 
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than 16 excisions), the system adapts by using the clutter cells to determine the overall 
system threshold level (V,,).. Thus, rather than excising the clutter, it is used to properly 
adjust the voltage threshold to maintain the false alarm rate. This is the same method 
emploved by the GO svstein. 

2. Extended Range Cells 

The hashed lines in Figure 55 show an additional 16 cells straddling the original 
16 leading and lagging reference cells. The additional cells are used only when a preset 
number of cells are excised from processing (indicating the possible presence of jammers) 
or when clutter edges dominate a reference window. This is the output of the binary 
integration. In general, only a marginal decrease in CFAR loss is obtained when going 
from a 32 to 64 cell system. This processing cost 1s well worth the effort however, when 
the system is under attack from multiple false target ECM svstems. Thus, when the false 
target jammung is detected or when immersed in clutter, a full 64 cell system is engaged 
into the EXGO processor to maintain robust performance. 

False target jamming (resulting in the 64 cell system) is declared by the excision 
logic once some predetermined number of excisions take place. For example if between 
eight and sixteen excisions take place among the original 32 leading and lagging cells, the 
system declares false target jamming. Since system performance of a standard 32 cell 
system is seriously degraded when high excision rates are used, the reference windows 
are expanded to the full 64 cell system. Thus, if every fourth cell contains a false target, 
the 64 cell system would excise 16 cells yet still maintain 48 cells for noise estimation. 

3. System Operation 

Other than the extended reference cells and the excision logic architecture, the 
EXGO processor behaves exactly asa GO CFAR. After the cells pass though the ex- 
cision logic, all remaining cells are summed by neighborhood and normalized by the 
proper number of non-interferer cells (nl, n2). The resulting values, yl and y2 are then 
input to the ‘GO’ logic for determination of the largest value. The ‘GO’ output is then 
multiplied by the threshold multiplier 7,, yielding a comparator threshold voltage, V,. 
The cell under test is then compared to this value. If the test cell’s amplitude exceeds 
jean a target is declared. 

When the system is in the clutter edge mode, the normal ’GO’ process contin- 
ues. Since the clutter cells are purposefully passed through the excision logic, they in- 
evitably become the ‘GO’ selection yielding the properly inflated adaptive threshold level 


that maintains the false alarm rate. 
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Figure 56. Excision Logic 


C. PERFORMANCE ANALYSIS 

As a member of the GO family, the EXGO detector maintains the key advantages 
of GO logic. The EXGO system has the advantage of low CFAR loss in the homoge- 
neous environment, but more importantly maintains a GO systems ability to control the 
false alarm rate in the presence of a clutter edge. In this analysis, we will show that the 
IEXGO system clearly overcomes the GO system’s inherent vulnerability to multiple 
target situations. The price paid for this significant improvement is a small increase in 
system complexity and an additional CFAR loss (~ 0.1 dB) in homogeneous environ- 
ments caused by excising legitimate noise samples. In addition, problems could occur 
when the processor is in a clutter region and multiple false targets appear. 

The input noise samples to the EXGO are normally distributed N(Q,1). The Monte 
Carlo probability of false alarm versus 7, is shown in Figure 57. Figure 58 shows the 
probability of detection curves using the threshold multipliers that offer false alarm rates 


of 107, 10-4, and 10%. Excision logic threshold multiplier (7,) was chosen to be three 
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times the running average noise level for these simulations This level led to an excision 
rate of approximately five percent of the input noise samples. 

Figure 59 shows the ability of the EXGO system to handle the clutter edge problem. 
Adaptive level 1 (excise level) and adaptive level 2 are detailed in this figure. As shown, 
adaptive level 1 1s higher than almost all the noise only peaks, vet is continually crossed 
in the clutter region. This graphically shows why excisions do not occur once the system 
recognizes that it 1s in a clutter edge. Adaptive level 2 shows the ability of the system 
to handle both the leading and lagging edges of this simulated sea clutter; thereby 
maintaining the false alarm rate. 

System operation in the multiple target environment 1s shown in Figures 60 to 65. 
In each of these plots, a P,, 10-* was chosen along with 7, = 3. All interfering targets 
have a SNR of 15 dB. In Figure 60, the EXGO svstem P, curves show relatively no 
change between the two interferer scenario and the noise only environment. Figure 61 
contrasts these curves by displaying the dramatic loss in capabilities of the standard GO 
CFAR with the same two interferers. For example, when the SNR of the primary target 
is 12 dB, the EXGO svstem maintains a detection probability of 0.72 whereas the GO 
CFAR yields a detection probability of 0.28. Figure 62 and 63 shows the loss in per- 
formance caused by four interfering targets and Figure 64 and 65 shows the performance 
caused by six interferers. In these plots the EXGO performance remains robust facing 
additional interferers whereas the GO system degrades significantly. 

To demonstrate the EXGO system in a multiple false target situation, 15 dB inter- 
feres Were injected into every fourth reference cel] window. In this situation, the ex- 
tended leading and lagging cells are engaged since between cight and sixteen excisions 
occur. This performance is shown in Figure 66. The resultant EXGO system has only 


and additional 0.2 dB CFAR loss whereas the GO CFAR performance 1s extremely poor. 


109 


PROBABILTY OF FALSE ALARM - EXGO CFAR 
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Figure 57. EXGO Probability of False Alarm Curve 
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EXGO CFAR PERFORMANCE IN CLUTTER EDGING 
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Figure 59. EXGO Performance in Clutter Edges 
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Pd Curves for EXGO CFAR -> Zero vs Two Interfering Targets 
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Figure 60. Effect of Two Interferers on EAGO CFAR 
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Pd Curves for EXGO CFAR -> Zero vs Four Interfering Targets 
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Figure 62. Effect of Four Interferers on EXGO CFAR 
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GO CFAR vs EXGO CFAR With Four Interfering Targets 
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igure 63. EXGO vs GO with Four Interfering Targets 
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Pd Curves for EXGO CFAR -> Zero vs Six Interfering Targets 
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Figure 64. Effect of Six Interferers on EXGO CFAR 
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GO CFAR vs EXGO CFAR With Six Interfering Targets 
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Figure 65. EXGO vs GO with Six Interfering Targets 
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EFFECTS OF FALSE TARGET JAMMING IN EXGO VS GO CFAR SYSTEMS 
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Figure 66. EXGO ys GO in Multiple False Target Jamming 
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D. SUMMARY 

The EXGO CFAR processor shows significant improvements over the conventional 
GO system. The price paid for these improvements are in the form of a small additional 
CFAR loss in the homogeneous environment due to excising a small number of legiti- | 
mate noise samples. Also, additional detector complexity is required as compared to the 
standard mean level detectors. Although complex, the EXGO system is easily imple- 
mented and performs faster than many of the rank ordering system that require cell 
sorting routines. 

The introduction of the extended reference cell concept enables a system to be 
adaptive to its real time operating situation. The ability to shift between 32 or 64 total 
reference cells enables the system to conserve its resources in the benign environment 
and increase its capability in large multiple target/false target jamming situations. As 
shown, even at relatively high excision rates, The EXGO adaptive performance is supe- 


r1or. 
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APPENDIX 


The four following programs are the GO CPAR and EXGO CIAR probability of 
detection and probability of false alarm curve generating code. In all programs an En- 


velope Approximation detector 1s used. 


12] 


[HEH HEHHRHHHHHHHHHOEDHHESEEEEEEHEEEEESCEEEEHEDESEDEDE DEED / 


/* THIS PROGRAM IS A MONTE CARLO SIMULATION OF A ’GO’ */ 


/* CFAR ALGORITHM. THIS PROGRAM PRODUCES THE PROB- ye 
/* ABILITY OF FALSE ALARM DATA POINTS TO BE MATLAB aes 
{> SE LOTIVED: ry 


[EOE E ELLIS SOP RENE A AE EO Oe See SP Se ee ere eee eae nana 
#include <stdio.h> 

#include <math.h> 

#include <stdlib.h> 

# nclude <time.h> 


[TEST HAH HHHAEHE THEE EH EEHEEHEDEHT HEHEHE HHEEEEDEETEE EE EEES / 


#fdefine PI (double) (4.0*atan(1.0)) 
#define TOTRV 1000000 


LEP OT ESTEE LES TELESIS SE NY SS MRS ASS A OSS SO SRS SO OLE ee eee, 


int z, first, overcount, supercount, randhold[2]; 

inka ,b yn; 

double threshold,a,b; 

double max, compval, rcdata, window([33],-numcellis; 

double percen,x1,x2,ul,u2, ldgval, ldgnorm, lagval, lagnorm; 


[eet ee sees ee eee eter eee eee See TORE OS TOS Se Oe Wee ee 8 28 88 88 8 8/, 


main () 
{ 
FILE. “writer; 


a=l; b=]; 
compval = 0.0; 
overcount=0; 
numcells=16.0; 
supercount=0; 


if ((writer=fopen("gopfa.dat”, "w") )==NULL) 
{ 
Print’ ne can de 7. 
exit (1); 
} 


LOT HOHHHH HHO HHERE TT IESE HEREET TT FERSEFESEREESESOSHEETERERED / 


/* THIS LOOP ITERATES THE THRESHOLD MULTIPLIER FROM 0 */ 
/* THOUGH 5.5. 100 DATA POINTS RESULT. wi 
LEPEPO SVP RS OSES PRESET PEACE CESAR OES POSS WOES O88 8 SOO 88 eee se / 
for (t=0;t<150; ++t) 

{ 

first=0; percen=0.0; 

overcount=0; supercount#=0; 

max=0 .0; 


threshold=((float)t)/10.0; 
srandom(1); 


[OP e PTA AAAS SEPA E TESS GSTS SECTS AOE COSTS OOS ES 8 Se eS eS ee, 


/* THIS LOOP ITERATES TRROUGHT THE 10 MILLION RANDOM 
/* VARIABLES TO GIVE ACCURACIES DOWN TO 10-7. 


[MARAE HSE HHH HHTEEE EEE EHET EEE HERE SETHE EH EEEHED ETE TED EDES / 


while (overcount<TOTRV) 
{ 
ldgval=0.0; lagval=0.0; 
for (i=0; 1<2; +41) 
{ 
randhold[i] = random(); 
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} 


[OPTS HSS SHH SHEESH HTH RE EHHEEEEEEEHHEHEHEHEHEEEEEETERERESD / 


/* BUILDING THE RANDOM VARIABLES TO BE UNIFORM. ne A 
JEP TSPE SETHE TEH OHHH HETHEEEHEHETERTEEERHEEHETEHEEEEEH EEE ED / 
ul=( (double) randhold[0))/2147483647; 

u2= ( (double) randhold[1}) /2147483647; 

mmeagrt ((—-2)*log(ul))*cos (2*PI*uz2) ; 

x2=sqgrt ((-2) *log(ul) ) *sin(2*PI*uz2) ; 

rcdata = fabs (x1) *a+fabs (x2) *b; 


enemies oe 6 EAS EAA E ROTH T ERS EE TEE ETE TERT ESE ET / 


/* INITIAL REFERENCE WINDOW LOAD UP. J +7 
eras Se FSF SSS C ESSEC EET ES TO EH EEE SEH ETHOETETEREHEE ES / 
if (overcount <33) 

{ 

window[32-overcount] = rcdata; 

overcount +1; 

} 
else 


[ETE TH ER HE RHE AH EHETEHHHHHEHHHHHHEEEEHEEEHEHHHEHEteates / 


/* SLIDING THE REF WINDOW AND INPUTTING A NEW R.V. * 
LOTTE HEHEHE HHH HH HEHEHE AAHEHHEEEATHHTHEEHEHEEEHEDEHEEEERE EES / 
{ 
for (n=0;n<32;++n) 

{ 

window [32-n] =window[31-n]; 

} 
window [0]=rcdata; 
overcount+=1; 


[LTTE T EHH HE THHH HE EH EHHEDEEAHHEETHHHEETEREHEREHEEEEEEEEEEESE / 


/* SUMMING THE LEADING AND LAGGING REF CELLS THEN gy 
/* NORMALIZING BY THE NUMBER OF CELLS. */ 
Peete ee RENCE TET ET HE HEAT HE EAEE HEHE EO TECETERETETESS / 
for (n=0;n<16;++4n) 

{ 

ldaval+=window([n); 

} 
for (n=17;n<33; +4+n) 

{ 

laqval+=window([n]; 

} 


laqnorm=lagval/numcells; 
ldgnorm=ldagval/numcells; 


[THO HEHE HEHE HEHE EE HEH EHEHEHOEEEHHHEEEEEHEHEEHESEHEHE DS / 


/* FINDING THE GREATEST OF VALUE AND USING IT TO DET- */ 
/* ERMINE THE ADAPTIVE THRESHOLD VALUE. THE GO IS MULT-*/ 
/* IPLIED WITH THE SCALING FACTOR. IF TRIS VALUE IS af 
/* LESS THAN THE TEST CELL, A TARGET IS DECLARED. ia f 
Perens SSS SAT STA SST OEE EEE TTT ESTE HE EEE ETHER EEE E OSS / 
if (ldqnorm>lagnorm) 

{ 

max=ldqnorm; 


else 
{ 
max=lagnorm; 


} 


compval=threshold*max; 
if (window[16)>compval) 
{ 
supercountt+=1; 


} 
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} 
} 
percen= ( (double) supercount) / (TOTRV-33) ; 
fprintf (writer, "$f %11.10f\n", threshold, percen) ; 
} 
fclose (writer); 
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PPPS HS CAE SHHFEAD OTOH OATES HEREC HES EASE H EEE H SSH SEHEHE TENA EE EEE S / 


/* THIS PROGRAM IS A MONTE CARLO SIMULATION OF A ’GO’ CFAR 7) 
/* ALGORITHM. THIS PROGRAM PRODUCES THE PROBABILITY OF DETECTION*/ 
/* DATA POINTS TO BE MATLAB PLOTTED. tf 


Seen eeeeeeee ee oo 2S 6 SSO AA ORE EL EE SERRA AER AEE EAE LEDS STE AE OES SS 


#include <stdio.h> 
#include <math.h> 


#tdefine NCELL 20000 
#define SAMPLE (double) (1.0/200.00) 
#define PI (double) (4.0*atan(1.0)) 


ai se ttt ORS HO OSS S / 


#define TOTRV 10000 
#define THRESHOLD 3.65 
#define SNRFROM -10.0 
#define SNRTO 30.0 


[eter eteteeeaaeeaeae / 


int i, j,idelay; 

double ppl [NCELL), pp2 [NCELL), pp4(NCELL], pp8 [NCELL], pp16[NCELL) ; 
Bouolesx, ccc, cC0O, cclt, cclb, ddl, dd2, pl, p2, sum; 

double sumvar, summean, summ, noisepwr; - 
int overcount, randhold[3],n,i, supercount,'t; 

double snr,amp,numcells,a,b, ldgval, lagval, gonorm; 
double var,x1,x2,u1,u2,u3,1[33],Q[33], adaptive, percen; 
double phi,s,cut,go,threshprime; 
double maxl,max2,max3,minl,min2,min3; | ; 


main () 
{ 
Blue wricver: 


a=1.0; b=1.0; adaptive=0.0; numcells=16.0; 
overcount=0; 
supercount=0; sumvar=0.0; 


threshprime=((numcells-1.0) *TBRESHOLD) / (numcel1S~-THRESHOLD) ; 


[ee PST PPP SHA PAETEE SSSR ERA EEHEEEHEEPE RHEE EHEHEEHHEHE TEESE EH EE SES / 


for (1=0; 1<NCELL; ++1) [ett e tessa teeeeeeeeeeeeenes / 
{ /* PDF for test cell noise*/ 
x= (double) i*SAMPLE; [HH eee eee eeeeteteeeeeeeeee / 


cc0=—pow (x, 2.0)/(2.0*pow(b,2.0)); 

CElt=pow (a, 2.0) *pow(x, 2.0); 
cclb=2.0*pow (b, 2.0) * (pow(a, 2.0) +pow(b,2.0)); 

Gee=cco+(Cclt/cclb); 

ddl=a*x/ (b*sqrt (2.0) *sqrt (pow (a,2.0)+pow(b,2.0))); 

dd2=b*x/ (a*sqrt (2.0) *sqrt (pow (a, 2.0)+pow(b,2.0))); 

pl=#a*b*texp (ccc) *sqrt (PI/2.0)*erf (ddl) / (sqrt (pow(a, 2.0)+pow(b,2.0))); 
p2=a*b*exp (ccc) *sqrt (PI/2.0) *erf (dd2)/ (sqrt (pow(a,2.0)+pow(b,2.0))); 
poma=(4.0/ (2.0*PI*Aa*b) ) * (pl+p2); 

summ+= (ppl [1) *SAMPLE) ; 

summean+= (ppl [i] *SAMPLE*x) ; 

sumvar+= (ppl[i] *SAMPLE*x*x) ; 

} 
noisepwr=sumvar-—pow (summean, 2.0); 


renee 22 2 << 2s ee ee eke eee eee eee eee ee eee eee eee eee ee Net / 


if ((writer=fopen("gopd.dat”, "w"))==NULL) 
{ 

printf ("bad"); 

) 
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LH SE EHS HERS H ED EH HEH EEEEHEHE HAHAH EHEEEEHEEEOHREHEESEHEREGEEEESHEEHED ORE ED / 


/* THIS LOOP ITERATES THROUGH THE SNR RATIOS USED TO CREATE TARGETS*/ 
/* IN THE CELL UNDER TEST. WE INCREMENT SNR BY 0.02 GB STEPS TO 
/* ACHIEVE A SMOOTH PLOT OF SNR vs PROBABILITY OF DETECTION. 


[PPP C ete erases e sere sweetest ewe e et Oe ee ete ee Oe ee eee, 


S=SNPFROM; 

while (s<=SNRTO) 
{ 

overcount=0; 
percen=0.0; 
supercount=0; 
go=0.0; 


[PPPS eeRs te ses ees ee sees ese cess Oe 888 eee eee eS ee ee 8 eee ee eee er 


/* RERE WE FIND TRE VALUE OF SOME TARGETS AMPLITUDE. 
/* EQUAL TO THE SQUARE ROOT OF TWO TIMES THE SNR TIMES TRE EST- 
/* IMATE OF THE NOISE POWER: THIS VALUE IS THEN ADDED TO THE 
/* NOISE ALREADY FOUND IN THE CELL UNDER TEST. 


| aT iD I SD Nia DE IAS TT EI ONE I 


noisepwr=1 .0; 
snr=pow (10.0, (8/10.0)); o 
amp=sqrt (2*noisepwr*snr); 


srandom(1); 5 we 


[PPS HPAAHIMH SSS IA EASELS SOOT S ESR O SET OS TRE SOOO 8 8 58 OD tS oe Oe See ee eee 


/* THIS LOOP ITERATES THROUGB ALL THE RANDOM VARIABLES. 


LETTS HFIP HESHENEF ONTO TH VOTH ODVERONIFETEFVE PEON CED OCTET EV OEE TEV EAS OOS O87) 


while (overcount<TOTRV) 
{ 
ldaval=0.0; lagqval=0.0; 
for (170; 1<3; ++1) 
{ 
randhold[i]=random(); 
} 


AMPLITUDE 1S*7 


LTE SAE SPOS SOP PEPPER SELES EE REP SS OES SSS SO 6 ORO RRR e Pe Oe eee Oe ens © Se ee eee 


/* BUILDING UP THE UNIFORM RANDOM VARIABLES 


ys 


ft EES POTOSI EES ES ESF Ee SO SOS 6 Oe Re OTe Re Oe eee ne nee eee eae 


ul=((double) randhoid[0))/2147483647; 
u2=((doubie) randhoid[{i])) /2147483647; 
u3=( (double) randhold[2)) /2147483647; 
x?=sart.({(—2) *10eq(u)))) “cos (2° Pi 712); 
x2=sqrt ((-2)* leg (1) ) *sin(2°Pi*u2); 


phi=(P1/4) *u3; 


FO a I EE ET SE ES TS DN SS I 


/* INITIAL REFERENCE WINDOW LOAD UP AND SLIDING OF THE REF CELLS. 


wae 


LAHORE PRET EREEHHEHR TED EEEEHEHEEERHEEHEHEOHEHEEEEERHEREEEEEDEEHESCHEDEEDED / 


1f (overcount<33) 
{ 
I {32-overcount ]=(x1); 
Q[32-overcount ] = (x2); 
overcountt+=1; 
} 
else 
{ 
for (n=0;n<32; ++n) 
{ 
I[32-nJ=I[31l-n]; 
Q(32-nJ=Q[31-n]); 
} 


IT{[O]=(x1); 
Q(0)=(x2); 
overcount+=1; 


Peer eee = oe ee eS ee OTP S PERO CCS H AER HSH H Re Ree nw ee eeweneewereree eer Tee / 


/* USING AN ENVELOPE DETECTOR, THE CELL UNDER TEST VALUE IS FOUND*/ 
eres 2 en ee eee SO eC ee eee eee eee eee eet eer ee eee eee eee eee eer es eee / 
if (fabs (I1(16}+ (amp*cos (phi) ))>fabs (Q([16}+ (amp*sin (phi) ))) 

{ 

maxl=fabs (1 (16}+(amp*cos(phi))); 

minl=fabs (Q(16)+(amp*sin(phi))); 

} 
else 

{ 

maxl=fabs (Q[16}+(amp*sin (phi) )); 

minl=fabs(1[16)+(amp*cos(phi))); 

} 


cut=a*maxl+b*minl; 


eeeneeeeeeeeerens = 2 = 2 2 Se SS AANA SSE T PES EH TOES E HEATH ERARETHEOAEAETEREE OHS ES / 


/* SUMMING THE LEADING AND LAGGING CELLS FOR ’GO’ DETERMINATION. yf 


nee we Oe ewe wee ee eee ee eevee eee eres seers ewer ere teeter eee / 


for (n=0;n<16;++n) 
{ 

if (fabs (I[n}))>fabs(Q{n))) 
{ 
max2=fabs(I[n)); 
min2=fabs (Q[n)); 
} 

else 
{ 
max2=fabs(Q{n)); 
min2=fabs(I[n}); 


ldoval+=a*max2+b*min2; 


} 


for (n=17;n<33; ++n) 
{ 
if (fabs (I[n})>fabs(Q[n})) 
{ 
max3=fabs (I[n}); 
min3=fabs(Q(n)); 
) 
else 
{ 
max3=fabs(Q[n)); 
min3=fabs(I{n)); 
} 
lagval+=a*max3+b*min3; 
} 


if (ldgqval>laqval) 
{ 

aqo=ldgval; 

} 
else 

{ 

go=lagval; 

} 


ee ee SSeS SASS HTH ATA E THES HEEAEATHE TRO EHECEEAEAEHEAREHESEPEETEREEE EEE ED / 


/* FINDING THE ADAPTIVE THRESHOLD LEVEL. IT IS THE PRODUCT OF THE i 
/* NORMALIZE ’GO’ AND THE SCALING FACTOR. THE CELL UNDER TEST IS THEN */ 
7~ COMPARED TO THIS LEVEL TO DETERMINE IF A TARGET IS DETECTED. yf 


eee e S22 2 SS 2 SSE R OKA E AAT ARETE OES OTE SOTTO N STEREO E WS TET OEE O RES / 


gonorm=go/numcells; 
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adept ive=THRESHOLD*gonorm; 


if (cut>=adaptive) 
{ 
Supercount+=1; 
} 
} 
} 
percen=( (double) supercount) / (TOTRV-31) ; 
fprintf(writer, "tf %11.10f£\n",s,percen) ; 
S=S+.2; 
} 
fclose(writer); 
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LA eee Pst SSGCTIREERSEHPEHERHHDEDRE HES HES ADORE EHTE REDS TE ERESE SEES Ee / 


/* THIS PROGRAM IS A MONTE CARLO SIMULATION OF A CFAR ALGORITHM */ 
/* THE ALGORITHM IS THE AUTHORS OWN AND IS CALLED ’EXGO’ WHICH */ 
/* IS SHORT FOR EXCISION GREATEST OF. THIS PROGRAM PRODUCES THE*/ 
/* PROBABILITY OF FALSE ALARM DATA POINTS TO BE MATLAB PLOTTED */ 
[POT HH HHH HEHE HATE DHE REE EEHE HEHEHE EHEHEAENHTEHEEH EEE SD OHEEHAHEH ED / 
#include <stdio.h> 

#include <math.h> 


#define EXCISELVL 3.0 
#define TOTALRV 1000000 
#define PI (double) (4.0*atan(1.0)) 


PET TAFT E FO CHTH TH AERAE EEE EEEEHED EFAS HEHEEEHEHEEEROEHOEEEEHHEDES / 


int overcount, supercount, curcount, ldgcount, lagcount, randhold{[2]; 
tec, a,Nn, ¥Y-Cluttcount; 


Gouble a,b,x1,x2,ul,u2,9g0,1[65}),Q[65), adaptive, threshmult, lagval; 
double ldgval,cut,curest, curestnorm, exciseval, Macey CIR e rN 
double percen,cluttsum, cluttlevel; 


daca c eden oc esststcerecsceeheotrevereesorurneorresesoer/ 


main () 
{ 
FILE *writer; 
a=1.0; b=1.0; 


if ((writer=fopen("far2.dat", "“w"))==NULL) 
{ 
printf ("fu") ; 
exit(1); 
} 


FTP TH HEH TATHHEHEDERHEEEHEEEUEDERTH HEP HERE E HEHEHE EHEHDHEDHEHEEEEEEEOEEHS / 


/* THIS LOOP ITERATES THE THRESHOLD MULTIPLIER FROM 0 THROUGE 5.5.*/ 
/* 550 DATA POINTS RESULT. 27 
oer eo CLS SSeS FEC CE EC EEP EERE OPO APE HEHE CASAC ECHPOPER EET OSH ERS  / 
for(t=0;t<126;++t) /* Threshold multiplier loop from T = 0 to 5.2 */ 

{ 

percen=0.0; 

overcount=0; 

curcount#=0; 

Ssupercount=0; 

go=0.0; 

y=0; 

curest=0.0; 

exciseval=10.0; 


threshmult=((float)t/25.0); 


srandom(1); 

[AAT HHHH HOHE OHHH HEHESERETHOHHOEHOHEEEH HEHEHE EEE ED EHHEEE OH EE EEE / 
/* THIS LOOP ITERATES THROUGH THE TEN MILLION RANDOM VARIABLES TO */ 
/* PRODUCE PFA DATA WITH AN ACCURACY DOWN TO 10-7. Lie d 


[THK HREM HHH EEA EHEE EHH EHHHEHEE EEE HEHEHE HEHE HHH HEHEHE EH EHH ERE EE EEE HES / 


while (overcount <TOTALRV) 
{ 
ldgval=0.0; 
laqval=0.0; 
lagcount=0; 
ldgcount=0; 


for (i=0; 1<2; ++i) 


{ 
randhold[i}=random(); 
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} 


‘ 
[OVC P SCSI ECHR SE CHESSER CRT EEE TS Oe Te Fe Fe eee Se Oe 


/* BUILDING THE RANDOM VARIABLES TO BE UNIFORM. = / 
fT ePe CREP EOE HPO eRe Tl Se 28 ere ee eee ee Se Oe ee Se 
ul=((double) randhold[0}) /2147483647; 

u2= ( (double) randhold[1)) /2147483647; 
xl=sqrt((—2)* log (ul) )*cos(2°*P1*u2) ; 

x2=sqrt ((-2) *log{ul))*sin(2*Fi*u2); 


[TETHER HHRHHHHHHHHHES HEHE HERE HEHEHE EE HOHE HE HEHEHE HEHEHE HEHEHE EDD / 


/* INITIAL REFERENCE WINDOW LOAD UP. */ 
fete tee ee eat OS Se Oe eee ee ee eee ee eee ere eee ieee een 
if (overcount<65) 
{ 
I{[64-overcount ]=xl; 
Q[64-overcount ) =x2; 
overcount+=]; 
) 
else 
{ 
YI li alae A ok oO Le Bi AN OL SO ES EN i th Ns A AN SL So IL ESD I IS DSI SIS IRIS CO 
/*SLIDING THE REFERENCE WINDOW AND INPUTTING THE NEW RANDOM 
/*VARIABLE. THE IF STATEMENT IS TRUE ONLY THROUGH THE FIRST 
/*REF WINDOW SLIDE. TRIS STATEMENT ALLOWS FOR THE INITIAL 
/*ESTIMATION OF THE CURRENT NOISE. ONCE INITIALIZED, THE EST- 
/*IMATE IS UPDATED WHEN A NEW RANDOM VARIABLE IS INTRODUCED 
/*INTO THE REFERENCE CELLS. 
LPH HHEHHHHHEHE HEHEHE EE EEHHEEEHEHEHEEEHHEEHEREEEH EHEC ED ERED EES 
for (n=0;n<64; ++n) 
{ 
T(64-nJ=#I[63-n); 
Q164-n)=#O[63-n); 
if (y==0) 
{ 
curest+=fabs (1[64-n)})+fabs(Q[64-n)); 
curcount +=1; 
} 
} 
y=1; 
T(Ol=xi> 
Q[Oj=x2; 
overcount+#=1; 


[AHHH HEHE HEHE HERE HEE HEHEHE HE DEH E HEHEHE TE ES EEEH EERE DERE EEE EE DS 


27 
es 
gi 
= 
i 
li & 
eed 
sat 


/ 


/* UPDATING THE NOISE ESTIMATE WITH THE NEW RANDOM VARIABLE say f 


/* ONLY'IF IS NOT RECOGNIZED AS AN ADDITIONAL TARGET OR JX. . 
[TPH Heat HHeHte et etee eee eee EHH HEHEHE RHEE EEE REE EEE EEE EE ED ED 
if (fabs (I1[0})+fabs (Q{0})<exciseval) 

{ 

curest+efabs(1[0})+fabs(Q[0}); 

curcount+=1; 

} 


[HHH HEHHHHHEH ED HEHEHE EHEHHHHHHS HHH EE EH HEED EHER EEE HEED EDD / 


/* EXCISEVAL IS DETERMINED TO BE THAT ADAPTIVE LEVEL IN */ 
/* WHICH THE REFERENCE CELLS ARE MEASURED AGAINST TO SEE */ 


/* IF THEY ARE NOISE VALUES OR POSSIBLE INTERFERING st 
/* TARGETS. EXCISE VALUE DETERMINED FROM THE CURRENT aA 
/* ESTIMATE OF TRE NORMALIZED NOISE. ard 


LTH HHHH HEHEHE RH HEHEHE HEHEHE HEHHEEEHTEEEEEEEETEFEENEH OEE TEDED / 


curestnorm=curest/curcount; 
exciseval=(curestnorm*EXCISELVL) ; 


[LOOT HHHHH HAHAHA HEHE EHO EEE HEHEHE REEEEEE EPH ETHER EHE DEDEDE TE / 


/* HERE THE LEADING AND LAGGING CELLS ARE SUMMED AS LONG */ 
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/ 
i 


/* AS THEY ARE LESS THAN THE EXCISE VALUE -ELSE THEY ARE */ 
/* IGNORED. =/ 


[OTP SO ES HH OHS HH HHH SEH EEHEEEHEEHEHEE SHEE OEHEREHEHOHHEDH EDO LED / 


Pon (n=) 6;n<32; ++n) 

{ 

if (exciseval>=fabs (I[n))+fabs(Q[n]l)) 
{ 
ldgval+=fabs (I(n))+fabs (Q(n))j; 
magcount +=] ; 
) 

} 


fOr (n=23;n<49; ++n) 

{ 

if (exciseval>=fabs(I[n])+fabs(Q[n))) 
{ 
lagqval+=fabs (I[n))+fabs(Q[n])); 
lagcount+=1,; 
} 

) 


[TPH HE HEE HE HERE EEE EE EERE EEE EEE EEE HEE EEE EEE ETHER E REESE HES / 


foe ize cooe checks for false target generating jamming vs 
/* If more than some percentage of targets are excised is 
for eiseessummed that FTG jamming iS occuring. In this */ 
/* case we assume that if 25% or > cellS are excised than */ 
/* we increase our excision threshold and expand our ref ‘*/ 
/* cell summation to 32632 in order to maintain ovr system*/ 
/* performance at high excision rates a 
[ETH HH HEHEHE HEHHEEHEEEHEEEEEEEEHEHEHEEEHE HE REEEHE HEHEHE BE REE / 
if(lagcount+ldgcount <= 24) 

{ 

lagcount=0; ldgcount=0; 

lagval=0.0; ldgval=0.0; 


if (lagcount+ldgcount <= 16) 
{ 
for (n=0; n< 64; +4+n) 
{ 
if(exciseval <= fabs(I[n})+fabs(Q[n})) 
{ 
eluttcount +#1; 
cluttsum+=fabs (I[n]))+fabs(Q[n]); 
} 

‘) 
cluttlevelecluttsum/ ( (double) cluttcount) ; 
exciseval=EXCISELVL*cluttlevel; 
cluttcount=0; 
cluttsum=0.0; 

) 
for (n=0;n<32;++n) 
{ 
if (exciseval>=fabs (I{n]))+fabs(Q[n]}) ) 
{ 
ldgval+=fabs (I[n))+fabs(Q[n]); 
ldgcount+=l; 
| ) 
) 


| for (n=33;n<64; ++n) 
{ 
if (exciseval>=fabs (I[n))+fabs(Q[n))) 
{ 
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lagqval+=fabs (I{n})+fabs(Q(n)); 
lagcountt#1; 
} 


} 


[TPIT HOSP ES CL ee eee eS oe fF es OS fe Oe eS SE © oe ee eee eee 


/* TRE NORMALIZED VALUES OF EITHER TSE 16X16 OR 32X32 ARE*/ 
/* DETERMINED. AT THIS POINT THE GREATER OF VALUE WILL he 
/* BE DETERMINED AND USED TO FIND TRE SYSTEM ADAPTIVE moh 
/* THRESHOLD LEVEL. aif 
[HHH HHH AH HE HE HHT TENTH HH HEHEHE HEHEHE HEHEH EEE ED ED EH EE HES / 
lagnorm=laqval/ ( (double) lagcount) ; 

ldgnorm=ldgval/((double) ldgcount); 


if (lagnorm>ldgnorm) 
{ 
go=lagnorm; 
} 
else 
{ - 
go=l]danorm; 
} 


adaptive=threshmult*go; 


[HHH HHH HHH HHEHETAHE SHE HHEEEHEDEHEDH HEHEHE HEHEHE ODEO DED EE EEE / 


/* KERE WE FIND THE VALUE OF THE CELL UNDER TEST AND COMPARE IT*/ 
/* TO THE ADAPTIVE THRESHOLD LEVEL. IF THE CUT IS GREATER WE */ 
/* INCREMENT SUPERCOUNT REFLECTING TARGET DETECTION. it 


[TTPO CAPES TES TEREST SPE PEL EDEVE CDSE TELS SS OS TSO SONS FS 8 Oe eile Or, 


cut=fabs (I1[32))+fabs (Q[32)); 


if (cut>adaptive) 
{ 
supercount+1; 
} 
} 
} 
percen=( (double) supercount ) / (TOTALPV-65) ; 
fprintf (writer, "$f %11.10f\n",threshmult, percen) ; 


} 
fclose (writer); 
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renee ee eat ee ees eS SSS SR Ee SS Se eS ee ee eee ee ee ee ee eer eee tere retro e ors eese re ee / 


Fs ck cach cc ce ee ca cs ce ee cee me ee a cr cre co re ee coe ee cee ee cr ee ce che elem ems ches rw com ae ee me me em ae ee hes ie i i ee oe i 
/* 

ioe XXXXXXXXXX XxX xx XXXXXXXXX XXXXXXXAXKXX 
ind xx XX XX xX xx xx 
‘ied XX xXx XX XX XX XX 
ed xx XX XX XX xX XX 
i XXXXXXXXX XXX XX XXXX xx XX 
Thad XXXXXXXXX x xX XXX XxX xX 
| oe xX XXX XX XX xX XX 
(oe xx xX" ex (7 xx x XX XX 
Pi XX XX XX : xx x XX XX 
Higa xX XX xX XX x XX XX 
| hig XXXXXXXXX XX XX XXXXXXXXXX XXXXXXXXAXX 
ae 
a a ee ee ee 


Pe SPS SSeS ES ASH SCHAAAASAFOCHENSHORSHHHESENSOEEDSREC DORSET ENE REO TED ETO EES / 


finclude <stdio.h> : 
#include <math.h> 


#define NCELL 20000 

#define SAMPLE (double) (1.0/200.00) 
#define PI (double) (4.0*atan(1.0)) 
eww er er erewrenses / 

#define TOTRV 5000 

#define THRESHOLD 2.42 

#define SNRFROM -10.0 

#define SNRTO 30.0 

#define EXCISELVL 3.0 


feeerteoearerawaaeone / 


[terete eeeteaetedredeeeeteeHetereeeeeeeDy / 

int i, j,idelay, bool, boolcount, cluttcount; 

double ppl (NCELL), pp2 (NCELL), pp4 [NCELL], pp8 [NCELL], pp16[NCELL] ; 
Hamee x, CCC, CCO, CC1t, cclb, ddl, GdZ, pl, p2, sum; 

Gouble sumvar, Summean, sumn,noisepwr, cluttsum, Cluttlevel; 


int overcount, randhold[3],n,supercount,t, curcount, ldgcount, lagcount, y; 
double snr,amp,a,b, ldgval, lagval, gonorm; 

Gouble xl,x2,ul,u2,u3,1(65),Q[65], adaptive, percen; 

double phi,s,cut,qo, curest, curestnorm, ldgnorm, lagnorm, exciseval; 
double maxl,max2,max3,minl,min2,min3; 


LTP HH HHH HEHEHE HRHHEEEERERHE ERE EHE DHE EES / 


main () 
{ 
FILE *writer; 


eee; b=] .0; 
overcount#0; 
supercount=0; 
bool=0; 

boolcount#=0; 


[TTS H HH HH HH SHH HHAEEHEHHRHEHHHEHREEHEEHH EE HHE REE EEE HEE RHEE E HEHEHE HEED / 


1f ((writer=fopen ("exgo2pda.dat", "w") }==NULL) 
{ 

peintt ("bad"); 

} 
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[PHOT SSOHHEHE DE SE HHO HEESHADEOCHHETETHEREEESESEHECEHTS OTERO EEWCHEDES © / 


ye 


THIS OUTER LOOP ITERATES THROUGH THE SNR VALUE RANGE SELECTED. */ 


[OPES OE HED OS Oe 8 tee eee ee ee re ee 


s=SNRF ROM; 
while (s<=SNRTO) 
{ 

overcount=0; 
percen=0.0; 
supercount=0; 
go=0.0; 

y=0; 
curcount=0; 
curest#0.0; 
exciseval=10.0; 


noisepwr=1.0; 
snr=pow(10.0, (s/10.0)); 
amp=sqrt (2*noisepwr*snr); 


srandom(1); 


[THEE AHHH HHH TEHEEH HEHEHE EHEEH ESOT EHTE HEE HEHEHE EHEEE EEE HES / 
e 


/* THIS INNER LOOP ITERATES THROUGH THE RANDOM VARIABLES * 7 
LPS EPEAT ER DEFOR ARTE EECA OPEEE AEE ESOC ERP SE SOROS ERS OF > 8 See ee ee 
while (overcount<TOTRV) 
{ 
ldgval=0.0; laqval=0.0; 
ldgcount=0; lagcount=0; 
for (i=0;1<3;++1) 
{ 
randhold[ijJ=random(); 
} 
ul=( (double) randhold[0})) /2147483647; 
u2=( (double) randhold[1}) /2147483647; 
u3=( (double) randhold[2)]) /2147483647; 
xl=sqrt ((-2) *log(ul)) *cos(2*PI*u2); 
x2=Sqrt((—2)*4100 (ul)) “sint2*PivuZ) 


phi=(PI/4)*u3; 


[Ot Ter ere eee eta e ae een tee eee 68 00 ee 6 ee 88 68S OF eee eee 


/* INITIAL REFERENCE WINDOW LOAD UP AND SLIDING ACTION sf 
[TTP HEHEHE HHHTERHEHDHOHOENEEOEETEEHEEEENEEEOCHEHEOEHEEEEEEEEESD / 
if (overcount<65) 

{ 

I [(64-overcount ]=(x1); 

Q[64-overcount )=(x2); 

overcountt+#1; 

} 
else 

{ 


for (n=0;n<64; ++n) 
{ 
I(64-n)=I[63-n]; 
Q(64-n)=Q[63-n); 
if (y==0) 
{ 
curest+=fabs (I [64-n))+fabs (Q[64-n])); 
eEnuy count +=1; 
} 
} 
y=1; 
I({O)=x1; 
Q(0)=x2; 
overcount+=l1; 
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JPET OHRHHER EH OHO EHH EHERERHEREHHETHHHEEEROHHEEEHHEH ENO ED & / 


/* USING THE NEW RV TO ADJUST THE EXCISE VALUE AVERAGE. */ 
SPECS HEEOEH EHH THR EHEHIREREEHEHEHEEHEHHEREHEAHEES EDEN E EE / 
if (fabs (1[0])+fabs(Q[0])<exciseval) 

{ 

curest+=fabs (1[0))+fabs(Q[0]}); 

curcount+=1; 

} 


ee SSP OH SESSA KE KEARSE REEA EEA KAA ERATE RRA EE HRS E e / 


/* NORMALIZING THE RUNNING NOISE ESTIMATE AND THEN i 
/* FINDING THE EXCISE VALUE BASED ON THAT ESTIMATE. as 


[THA HHHHEHHEHHHHHHHHEEE HEHEHE EEHE TENE EERE ERE EREEEED EES / 


curestnorm=curest/curcount; 
exciseval=(curestnorm*EXCISELVL) ; 


JAAP EHH HHH HHH HHAEEEEAEEHEETEAEEEEEWEEEEEEEEEHHEEEHEHE RHE EEE SH / 


/* summing the lagging and leading ref cells that do not */ 
/* exceed the excision value. Bool is used to count the */ 
feeeeeal number is excisions for the first run. This data*/ 
/* will help us decide on a proper EXCISELVL to choose 80 */ 
/* that very few data points are excised in a normal noise*/ 


/* environment +] 
a Se SS ee SS AE SSAA TPT SASSER AEA HEHEHE TE HSER E REET ES / 


for (n=16;n<32; ++n) 
ce aion casi ini) <eaneiaih 
Pee sehr inj) txts Orn); 
faecount+=1; 
ee 
eee sic0) 
bool count +=4; 


} 


for (n=33;n<49; ++n) 
{ 
if (exciseval>=fabs(I[n])+fabs(Q[n]})) 
{ 
lagval+=fabs (I[n]))+fabs (Q[n}); 
lagcount+=1; 
} 
else 
{ 
if (bool==0) 
{ 
boolcount+=1; 
} 


} 


{ttt Hee aeeeneeeeeeeene ee eeR EEE HEH HEE HEHEHE EH EEE ERE EEE EE Ee Ee / 


/* False Target Jamming is checked here. If greater than a / 
/* 25% of the data samples are excise FTG jamming is assumed */ 
/* We then reset the ldg and lagging sums to zero and recalc-*/ 
/* ulate based on 32X32 ref cells and a higher excision level*/ 


a SX SRS SAATE SERS E ESE RHE SERS ASSERTS AA HE SEA SSA eA See eee Se / 
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if (lagqcount+ldgcount <#=24) 
{ 
lagcount=0; ldgcount <0; 
ldgval=-0.0; lagval=0.0; 


if (lagcount +ldgcount <<16) 
Se eer eer 
if (oicl se yalcceeecen a eee 
cece 
cluttsum+<fabs (I[n))+fabs(Q{[n}); 
' } 


cluttlevelecluttsum/ ( (double) cluttcount) ; 
exciseval=EXCISELVL’cluttlevel; 
cluttcount<0; cluttsum=<0.0; 

} 


for (n=0;n<32; ++n) 
{ 
if (exciseval>«fabs (I[{[n})+fabs(Q[n))) 
ldgval+=fabs (I[{n))+fabs(Q([(n}); 
ldgcount+#1; 
) 
} 
for (n=33;n<65; +4+n) 
{ 
if (exciseval>=fabs (I({[n})+fabs(Q{n))) 
{ 
lagval+«fabs (I[{n))+fabs(Q[n)); 
lagcount+#1; 
} 
} 
) 


[LOPS HHH OTH EHEEEHEOHE OLEH EEHEEEEEEHEOEEREHEEREEEEEHEEEEEDSE / 


/* NORMALIZING THE LAG AND LEADING REF WINDOWS THEN as 
/* CHOOSING THE GREATER OF. THIS VALUE IS THEN USED TO */ 
/* DETERMINE THE ADAPTIVE THRESHOLD LEVEL. */ 


[FSSA SHEER HHHEHHHOEEOHOEREEHEHEEEREE HEHEHE EHEEEEEEHEHE ETE EE / 


lagnorm=lagval/ ((double) lagcount) ; 
ldgnorm=ldgval/ ( (double) ldgcount) ; 


if (lagnorm>ldgnorm) 
ee 

ee 
Ate Rar 
} 


adapt ive=THRESHOLD*go; 


LEHRER EHERE HH OHH TEEHEEE HOE EE EETHEEHHHREEHEOHEEEEHEETED ESE / 


/* KRERE WE FIND THE VALUE OF TRE TEST CELL AND THEN COMP-*/ 
/* ARE IT TO THE ADAPTIVE THRESHOLD LEVEL FOR DETECTION. */ 


[PTH HAHOHH HAHAH H HH EET EH EEEEEHEHETHEHEHeEHHeNetteeetenetee / 


cut=fabs (I [32]+(amp*cos (phi) ))+fabs (Q{[32)+(amp*sin(phi))); 
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if (cut>adaptive) 
{ 
supercount+=1; 


} 
} 
} 
percen=((double) supercount) / (TOTRV-65) ; 
forint tf (writer, "%£ %11.10f£\n",s, percen) ; 
S=S+.2; 


bool+=1; 


} 
fclose (writer); 
perntt ("td", boolcount) ; 
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